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Abstract 



Our aim is to analyze the various energy functionals appearing in 
the physics hterature and describing the behavior of a Bose-Einstein 
condensate in an optical lattice. We want to justify the use of some 
reduced models. For that purpose, we will use the semi-classical anal- 
ysis developed for linear problems related to the Schrodinger operator 
with periodic potential or multiple wells potentials. We justify, in 
some asymptotic regimes, the reduction to low dimensional problems 
and analyze the reduced problems. 
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1 Introduction 



1.1 The physical motivation for Bose-Einstein conden- 
sates in optical lattices 

Superfluidity and superconductivity are two spectacular manifestations of 
quantum mechanics at the macroscopic scale. Among their striking charac- 
teristics is the existence of vortices with quantized circulation. The physics 
of such vortices is of tremendous importance in the field of quantum fluids 
and extends beyond condensed matter physics. The advantage of ultracold 
gaseous Bose-Einstein condensates is to allow tests in the laboratory to study 
various aspects of macroscopic quantum physics. There is a large body of 
research, both experimental, theoretical and mathematical on vortices in 
Bose-Einstein condensates |PeSm[ [PiSt[ [M \ iLSSYj . Current physical inter- 
est is in the investigation of very small atomic assemblies, for which one 
would have one vortex per particle, which is a challenge in terms of detec- 
tion and signal analysis. An appealing option consists in parallelizing the 
study, by producing simultaneously a large number of micro-BECs rotating 
at the various nodes of an optical lattice [Sn]. Experiments are under way. 
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A major topic is the transition from a Mott insulator phase to a superfluid 
phase. We refer to the paper of Zwerger [Z] and the references therein for 
more details. Our framework of study will be in the mean field regime where 
the condensate can be described by a Gross Pitaevskii type energy with a 
term modeling the optical lattice potential. The mean field description of 
a condensate by the Gross Pitaevskii energy has been derived as the limit 
of the hamiltonian for N bosons, when N tends to infinity [LSYt ILSj in the 
case of a condensate without optical lattice. The scattering length of the 
interaction in the iV-body problem is such that Na^ g. The rigorous 
derivation in the case of an optical lattice where there are fewer atoms per 
site is nevertheless open. In a one-dimensional optical lattice, the condensate 
splits into a stack of weakly-coupled disk-shaped condensates, which leads to 
some intriguing analogues with high-Tc superconductors due to their similar 
layered structure \SnStl\ [SlISt2l IKMPSl lABBli IABB21 IXBS] . Our aim, in 
this paper, is to address mathematical models that describe a BEG in an 
optical lattice. The theory which we will develop is inspired by a series of 



physics papers [Snl IS^IStTl IS^I5t2l IKMPSII5TKB] . We want to justify their 



reduction to simpler energy functionals in certain regimes of parameters and 
in particular understand the ground state energy. 

The ground state energy of a rotating Bose-Einstein condensate is given 
by the minimization of 




(1.1) 



under the constraint 




(1.2) 



where 



• = + , r = (x, y, z) , 



• r2 > is the rotational velocity along the z axis. 



• X r = Q{—y, X, 0) , 



• > is the scattering length. 
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The experimental device leading to the realization of optical lattices requires 
a trapping potential V{r) given by 

nr) = ^Kr^ + -.V), (1.3) 

corresponding to the magnetic trap. We assume that the radial trapping 
frequency is much larger than the axial trapping frequency, that is 

0<u,«uj±. (1.4) 

We will always assume the condition : 

0<n<uj± (1.5) 

for the existence of a minimizer. The trapping has to be stronger than the 
centrifugal force. The presence of the one dimensional optical lattice in the 
z direction is modeled by 

W,iz) = \w{z), (1.6) 

where \ is the lattice deptlfl, and w is a positive T-periodic function. In the 
whole paper, we will assume : 

Assumption 1.1. 

The potential w is a C°° even, non negative function on M which is T-periodic 
and admits as unique minima the points kT (k & Z). Moreover these minima 
are non degenerate. Thus, 

w{z + T)= w{z) , w(0) = , w"(0) > , w{z) > if z ^ TZ . (1.7) 



An example is 

w{z) = sm\^) (1.8) 
A 

and A is the wavelength of the laser light. The optical potential creates 
a one- dimensional lattice of wells separated by a distance T = A/2. We 
will assume that e tends to (this means deep lattice) and that T is fixed. 
Furthermore, we assume that the lattice is deep enough so that it dominates 
over the magnetic trapping potential in the z direction and that the number 
of sites is large. Thus we will, in this paper, ignore the magnetic trap in the 
z direction, and this will correspond to the case 

u, = 0. (1.9) 

^called Vz in Snoek [Sn] 
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We will mainly discuss, instead of a problem in M^, a periodic problem in the 
z direction, that is in x [— |^[, where T corresponds to the period of the 
optical lattice, or in M^,^ x [— ^[ for a fixed integer iV > 1. Therefore, 
we focus (see however Subsection 18.21 for a justification of this choice) on the 
minimization of the functional 



NT NT r 
" 2 ' 2 i 



iVt X r^l^ - -Vt'r^\mY + {V{y) + W,{z))\mY + ^|*| 



under the constraint 



/ 

Jm 



|\E'(x, -2)P dxdydz 



with 



nr) = 2^±V, (1.12) 

the potential given by fll.6p -( fTT7j) . and the wave function satisfying 

^{x,y,z + NT) = ^!{x,y,z) . (1.13) 

This functional has a minimizer in the unit sphere of its natural form domain 
(see (18.101) for its description) S^^'^ and we call 



N 



inf gr^'l^)- 



■per,N 



;i.i4) 



dxdydz , 
(1.10) 



Notation 

In the case = 1, we will write more simply 



E. 



,per,(N=l) 



When = 0, we will sometimes omit the reference to Q. 



;i.i5) 



Our aim is to justify that the ground state energy can be well approx- 
imated by the study of simpler models introduced in physics papers [Snl 
ISnStlt IKMPSj and to measure the error which is done in the approxiamtion. 
For that purpose, we will describe how, in certain regimes, the semi-classical 
analysis developed for linear problems related to the Schrodinger operator 
with periodic potential or multiple wells potentials is relevant: Outassourt 
|0u] . Helffer-Sjostrand [Hel [DiSj] or for an alternative approach jSi] . 
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1.2 The linear model 



The linear model which appears naturally is a selfadjoint realization associ- 
ated with the differential operator : 

Hn = Hl + H,, (1.16) 

with 

H'l := -^A,,, + ^ulr^ - ilL, , (1.17) 

= i{xdy - yd^) , (1.18) 

and 

H.:=-l-^ + W,iz). (1.19) 

In the transverse direction, we will consider the unique natural selfadjoint 
extension in L^(R^j^) of the positive operator by keeping the same no- 
tation. In the longitudinal direction, we will consider specific realizations of 
Hz and in particular the T-periodic problem or more generally the (NT)- 
periodic problem attached to which will be denoted by H^'^^ and H^'^^'^ 
and we keep the notation for the problem on the whole line. 
So our model will be the self-adjoint operator 

^per,N + RP/''^ . (1.20) 



In this situation with separate variables, we can split the spectral analysis, 
the spectrum of H^^'^ being the closed set 

a(i/r^) := a(i/^) + a(i/r'^) . (1.21) 



The first operator Hj^ is a harmonic oscillator with discrete spectrum 
as we will explain in Section O Under Condition (11.51) . the bottom of its 
spectrum is given by 

:= mi{a{Hl)) = , (1.22) 

hence is independent of Vt. 

A corresponding groundstate is the Gaussian = (^) ^ exp— ^r^. The 
gap between the ground state energy and the second eigenvalue (which has 
multiplicity 1 or 2) is given by 

5^ ■= \^^^- =uj^-n. (1.23) 
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The properties of the periodic Hamihonian H^'^^'^ , which will be de- 
scribed in Subsection 13.21 (Formulas (13.81) and (13. 9p for the physical model), 
depend on the value of A^. In the case = 1, we call the groundstate of if^^^ 
01 (2;) and the ground energy (or lowest eigenvalue) Ai,^. In the semi-classical 
regime e — > 0, Ai,^ satisfies 

Ai,. ~ ^, (1.24) 

for some c > 0. The splitting 6z between the groundstate energy and the 
first excited eigenvalue satisfies 

5. ~ ^ , (1.25) 

for some c > 0. 

For > 1, the groundstate energy of H^^^'^ is unchanged and the corre- 
sponding groundstate is the periodic extension of 0i considered as an 
(A^T)-periodic function. More precisely, in order to have the LF'- normaliza- 
tions, the relation is 

0f = ^0i, (1-26) 

on the line. But we have now exponentially close eigenvalues of the order 
of Ai_2 lying in the first band of the spectrum of Hz on the whole line. They 
are separated from the (A^ + l)-th by a splitting 5^ which satisfies : 

5f = 4 + a(exp -S/t)) . (1.27) 

Here the notation (9(exp —S/e)) means 

d(exp -S/e)) = C(exp -S'/e) , WS' < S . (1.28) 

The first A^ eigenfunctions satisfy 

cP^z + T)= eM^^^^^) 0f (^) , for £ = 1, . . . , AT , (1.29) 

corresponding to the special values k = '^'^^^^^ of what will be called later a 
/c-Floquet condition (see ( 1A.2P ). 

We will sometimes use another real orthonormal basis (called (NT)- 
periodic Wannier functions basis) {ip^) (j = 0, . . . , A^ — 1) of the spectral 
space attached to the first A^ eigenvalues. Each of these (A^T)-periodic func- 
tions have the advantage to be localized (as e — > 0) in a specific well of 
considered as defined on M/(A^T)Z. 
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1.3 The reduced functionals 

We want to prove the reduction to lower dimensional functionals by using 
the spectral analysis of the linear problem. There are two natural ideas to 
compute upper bounds, and thus find these functionals. We can 

• either use test functions of the type 



where ijj± is the first normalized eigenfunction of H^^ and minimize 
among all possible L^-normalized ^(z) to obtain a ID- longitudinal re- 
duced problem, 



where (f)i is the first eigenfunction of i/f^'" and minimize among all 
possible L^-normalized tp{x,y) to obtain a 2D transverse reduced 
problem, 

— or in the case A'^ > 1 



'^{x,y,z) = (f){z)ip±{x,y) , 



(1.30) 



• or use 



— in the case = 1, 



*(x,y,z) = (f)i{z)i/j{x,y) 



(1.31) 



N-l 




(1.32) 



j=0 



where il>^{z) is the orthonormal basis of Wannier functions men- 
tioned above, and minimize on the suitably normalized '0j,±'s 
which provide N coupled problems. Wc denote by IIjv the projec- 
tion on this space. For "if e L^{R^x] - ^[) , we have 



N-l 




(1.33) 



j=0 



with 
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Computing the energy of a test function of type fll.30p . we get 

gr'''W = ^± + ^A(0) (1-34) 

where is the functional on the A^T-periodic functions in the z direction, 
defined on H\R/NTZ) by 

NT , . 

with 

g ■■= g (^j Ji^±{x,y)\'^dxdy^ = ^guj±. (1.36) 

The functional is introduced by [KMPS] who analyze a particular case. 
Its study in the small e limit is one of the aims of this paper. 

For test functions of type fll.3ip . we get in the case N = 1 

Q^^^v]/) = Ai,, +£B,nW (1.37) 

with 



:i.38) 



and 

2 



g-=gyj ^\M^)rdz\ . (1.39) 

In the case > 1, we define ^^^n((V^i,-L)i=o,...,Af-i) by 

j 

with 

N-l 

vi> = 5^^f(z)^,-^(x,i/). (1.41) 
i=o 

Of course when minimizing over normalized ^P's, one gets more simply the 
problem of minimizing 

Qr''W = K^ + ^B,nm,±))- (1-42) 
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As such, the energy S^^ does not provide N coupled problems but one 
single energy depending on N test functions. Nevertheless, in the small e 
limit, the Wannier functions are localized in each well. Thus each function 
ipj^x only interacts with its nearest neighbors and this simplification pro- 
vides coupled problems, as suggested by Snoek |Sn] on the basis of formal 
computations. We will analyze their validity. This reduced functional is 
somehow related to the Lawrence-Doniach model for superconductors (see 
|ABB1[RBB2] ). 



1.4 Main results 

1.4.1 Universal estimates and applications 

The analysis of the linear case immediately leads to the following trivial 
and universal inequalities (which are valid for any A^ and any Q such that 
0<n<uj±) 

Xi,z + co±< ^r"^ < Ai,. + co^ + In, (1.43) 



where 



This universal estimate is obtained by using the test function 

v^P^^'^(x,i/,^) = V'±(x,l/)0f(^), 

where (p^ is the A^-th normalized ground state introduced in (I1.26P and 
4'±{x,y) is the ground state of H^, actually independent of Q. 
From fll.26p . we have : 

NT NT T 

" {cl^^[z)Ydz=^^ [ ' Mzrdz = ^ r Mz^dz, (1.45) 

NT l\l I NT I\ T 

where, as e — 0, and, under Assumption (11. 7p . it can be proved (see fl3.10p ) 
that 

T 
9 

4 



x{zf dz ^ Cit~-^ , (1.46) 



for some explicitly computable constant C4 > 0. 
Thus, we have the following estimate for J^v 



C4 gijJi. _i 

(1.47) 
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An immediate analysis shows that Ai^^ + uj±_ is a good asymptotic of 
^per,Af ^j^^ hmit e when g is sufficiently small (what we can call the 
quasi-linear situation). More precisely, we have 

Theorem 1.2. 

Under the condition that either 

(QLa) g , (1.48) 

or 

{QLh) guj^e^«l, (1.49) 

then we have 

i^r'^ = (Ai,. + a;^) (1 + 0(1)), (1.50) 

as e tends to 0. 

Each of these conditions implies indeed that /at is small relatively to 
or to oj^. 



So our goal is 

• to have more accurate estimates than (11.501) . 

• to analyze more interesting cases when none of these two conditions is 
satisfied and to give natural sufficient conditions allowing the analysis 
of reduced models. 

We are able to justify the reductions to the lower dimensional functionals 
£^ and E^o, when their infimum is much smaller than the gap between the 
ffist two excited states of the linear problem in the other direction, namely 
in case A, when is much smaller than where 

m^= inf £^(0), (1.51) 

and in case B, when is much smaller than 1/e, the gap between the two 
ffist bands of the periodic problem on the line, where 



m 



N _ i^t cN 

B,n 



inf,,^ S^M,,^)). (1.52) 



An independent difficulty is then to have more accurate estimates and 
m^Q according to the regime of parameters. We do not have universal esti- 
mates for this but have to separate two cases : 
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• the Weak Interaction case, where the interaction term (L"^ term) is at 
most of the same order as the ground state of the hnear problem in the 
same direction; 

• the Thomas Fermi case, where the kinetic energy term is much smaller 
than the potential and interaction terms. 

In what follows, when is not mentioned in m^, m'^^, S^, £^0.1 then 
the notations are for = 1. Similarly, if VL is not mentioned, this means 
that either the considered quantity is independent of VL or that we are treat- 
ing the case f2 = 0. To mention the dependence on other parameters, we 
will sometimes explictly write this dependence like for example m^(e,^) or 



1.4.2 Case {A) : the longitudinal model 

We consider states which are of type fll.30p with ip G L'^iM.z/ {NT)"!,). The 
energy of such test functions provides the upper bound 

£^r'''<^±+^A(e,?) (1.53) 

where is given by fll.511) and ^ was introduced in fll.36p . 

In order to show that the upper bound is an approximate lower bound, 
we first address the 'Weak Interaction" case, 

{AWIa) 1 << e{uj^ - n) , (1.54) 

and, for a given c > 0, 

{AWIh) giu±e'2<c. (1.55) 

The first assumption implies that the lowest eigenvalue Ai ^ of the linear 
problem in the z direction (having in mind fll.24p ) is much smaller than the 
gap in the transverse direction S± = uj± — Q. This will allow the projection 
onto the subspace ip± ^L'^i'^z/ {NT)Z). The second assumption implies that 
the nonlinear term (of order guj±/ ^/e) is of the same order as Xi^z- It implies 
using fll.24p . fll.47p and the universal estimate 

Ai,. < < Ai,, + , (1.56) 

that 




(1.57) 
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Here ~ means "of the same order" in the considered regime of parameters. 
More precisely we mean by writing fll.571) that, for any eg > 0, there exists 
C > such that, for all e G]0,eo], any g,uj_[_ satisfying (11.551) . 

Note that most of the time, we will not control the constant with respect to 
N. 

All these rough estimates are obtained by rather elementary semi-classical 
methods which are recalled in Section El More precise asymptotics of 
will be given under the additional Assumption fll.48p in Section 15. 2[ Thus, 
by (11.541) . is much smaller than 5±. We will prove 

Theorem 1.3. 

When e tends to 0, and under Conditions \l-54\) and M.55\) . we have 

Er'' = ^± + "^A(e,?) (1 + 0(1)). (1.58) 

We now describe the "Thomas-Fermi" regime, where we can also justify the 
reduction to the longitudinal model. We assume that, for some given c > , 

{ATFa) guj^^e»l, (1.59) 
{ATFh) guj±e'^<c, (1.60) 

(ATFc) g-^e-^uj^^^ « {uj^-n)l . (1.61) 

Note that ffT3^ is the converse of ffT33D while ffT3^ and ffTUTD imply that 
1 << e{LJ± — Q). This implies Ai^^ << 6±, which is the main condition to 
reduce to case A. Assumptions (I1.59P and (11.601) allow to show that : 

(1.62) 

and this also implies that the nonlinear term is much bigger than Sz- 

The estimate (11.621) will be shown in Section 15.31 together with more precise 

ones with stronger hypotheses (see Assumption (I5.19P and (I5.20p ). 

Theorem 1.4. 

When e tends to 0, and under Conditions (11.590 . (I1.60p and (ll.6ip . we have, 
as e —>■ 0, 

^r'' = ^±+<(e,?)(i+«(i))- (1-63) 

The proofs give actually much stronger results. 
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1.4.3 Case (B) : the transverse model 



This corresponds to the idea of a reduction on the ground eigenspace in the 
z variable, where the interaction term is kept in the transverse problem: 
therefore, this is a regime where u!±e << 1. We recall that we denote by Ai^^ 
the (iV-independent) ground state energy of H^^^'^ and by (p^ the normalized 
ground state. We consider states which are of type fll.3ip or (11.321) . We have 
defined S^^ by (ll.40p - (ll.4ip and m'^^, the infimum of the energy of such 
test functions by (11.520 . We have the upper bound 

Er"" < Am + <n ■ (1-64) 

When = 1, mB,Q is a function of g and uj± as it is clear from (I1.38P and 
(fL52l) . Note that, from ffLlHll . we get 

T 

~g = g{[' MzYdz)^^. (1.65) 

Again we can discuss two different cases according to the size of the interac- 
tion. In the Weak Interaction case, we prove the following : 

Theorem 1.5. 

When e tends to 0, and under the conditions 

(BWIa) ge-^2<C, (1.66) 

(BWIb) uj^e«l, (1.67) 

then 

i?r'' = AM + m^o(l + o(l)). (1.68) 



Condition (BWIb) implies that the bottom of the spectrum of the linear 
problem in the x — y direction is much smaller than 6z, the gap in the z 
direction, which is of order 1/e. Condition (11.660 . together with (11.430 and 
(ll.47p . implies that m^Q satisfies 

<,n ^ ^± ■ (1-69) 

Indeed, (BWIa) and (BWIb) imply ge^uj± « 1, that is (QLb). 
In the Thomas-Fermi case, we prove the following : 
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Theorem 1.6. 

When e tends to 0, and under the conditions 

(BTFa) ^fe«g, (1.70) 

{BTFh) u^^ei«l, (1.71) 

and 

(BTFc) gh-^uj^«l, (1.72) 

then 

Er'' = Ai,. + m^,^(l + o(l)). (1.73) 

Note that (BTFa) is the converse of (BWIa). We will see in Proposition 
16.31 (together with (16.41) . fl6.16p and fl6.49p ) that, under these assumptions 
and Assumption fl6.15p . the term m^Q satisfies 

m^^^^u^^/e'/\ (1.74) 
and thus is much smaller than 6^ which is of order -. 

Our proofs are made up of two parts : rough or accurate estimates of 
and m^Q^ on the one hand and a lower bound for E^^'^ on the other 
hand. The lower bound consists in showing that the upperbound obtained by 
projecting on the special states introduced above in (11.301) . (I1.3ip or (ll.32p 
is actually also asymptotically a good lower bound. 



1.4.4 Tunneling effect and discrete model 

Since the Wannier functions are localized in the z variable, the energy of 
a function ^ = J2j'=o''Pfi^)'4'j,±i^^y) provides at leading order the sum of 
decoupled energies for i/jj^± on each slice j. At the next order, in the 
computation of the norm of the gradient, only the nearest neighbors in z 
interact through an exponentially small term, describing what is called the 
tunneling effect. These simplifications are discussed in section [71 We are 
lead to new functionals and in particular a discrete model that we analyze 
in relationship with the physics papers. 

In case A, the behavior on each slice j is the same, given by ipx and it is 
the behavior on the z direction which has a tunneling contribution. There 
are no vortices whatever the velocity Q. 

In case B, for A^ = 1, there are vortices for large velocity and they are 
located on each slice at the same place. For A^ large, it is an open and 
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interesting question to analyze whether it is possible for a vortex line to vary 
location according to the slice, whether vortices interact between the slices 
and how. This could be performed using our reduced models. 

1.5 Organization of the paper 

The paper is organized as follows. In Section O we start the spectral analysis 
of the linear problems in the longitudinal and transverse directions. We recall 
in particular the main techniques which can be used for the analysis of the 
spectral problem with periodic potential on the line. Section [3] is devoted to 
the semi-classical results for the periodic problem. Although we are mainly 
interested in ID-problems we recall here techniques which are true in any 
dimension and can be useful for the analysis of 2D or 3D optical lattices. 

In Section HJ we prove the main theorems for case A. In Section [5l we 
analyze the ground state of the ID nonlinear energy £^ for = 1 and 
A^ > 1 and also distinguish between the two cases: Weak Interaction and 
Thomas- Fermi. Section [6] corresponds to a similar analysis for the transverse 
models £^ . Section [7] is devoted to the tunneling effects and discuss, on the 
basis of the semi-classical estimates of Section [3l some results by physicists 
on the discrete nonlinear Schrodinger model. In Section [HI we analyze various 
boundary conditions and compare in particular the problems on (which 
is completely solved) and the problems on M^x] — ^,^[ with periodic 
condition which seem physically more interesting. 

2 Analysis of the linear model 

The linear model which appears naturally is associated to 

Hn = + Hz , 

which was presented in the introduction (see f ll.17p - p.2ip ). A natural con- 
dition (for the strict positivity of the operator H^) is Condition (11. 5p . In 
this situation with separate variables, we can split the spectral analysis in 
the separate spectral analysis of and the spectral analysis of a suitable 
realization of Hz which will be presented in the next subsection. 
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2.1 The harmonic oscillator in the transverse variable 

For simplicity, we begin the analysis of H± = with the case = 0. 
The first operator is a harmonic oscillator with discrete spectrum and 

the bottom of its spectrum is given by 

mfiaiHl)) = u;^. (2.1) 

A corresponding L^-normalized ground state is the Gaussian 

^.= (^^)^xp-^r^ (2.2) 

Moreover the gap between the ground state energy and the second eigenvalue 
(which has multiplicity 2) is given by 

A2,± -uj±=u±. (2.3) 

The spectrum of can be recovered by considering first the joint spec- 
trum of and L^. For each eigenspace of Lz corresponding to i for some 
£ G Z, we can look at the operator 

Hf^ ■.= ~A,,y + ^u;^\'-ni, (2.4) 

considered as an unbounded operator on L'^(M?) fl Ker (L^ — i). 

More precisely, for Q satisfying 01. 5p . a common eigenbasis of and 
is given by the set of (not normalized) Hermite functions: 

,(x, y) = e'^(-'+^') (d, + idyY {d, - idyf {e-^^i-'+y')'^ (2.5) 

where j and k are non-negative integers. 
The eigenvalues are {j — k) for and 

Ej^k = cu± + - + {uja. + ^)k (2.6) 

for H^. 

The spectrum of H^^ is obtained by considering the pairs {j, k) such that 
j-k = L 

We emphasize that this orthogonal basis of eigenfunctions is independent of 

n. 
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2.2 The band spectrum in the longitudinal direction 



The second operator can be analyzed by semi-classical methods but note 
that our semi-classical parameter is e. One can of course in the case of the 
specific w introduced in (11.71) recognize this operator as the Mathieu operator 
(for which a lot of information can be obtained using special functions (see 
|ASj )) but we prefer to give the presentation of the theory for a more general 
periodic potential w. We hope that the general ideas which are behind will 
become clearer. 

There are two related approaches for the analysis of the spectrum of 
Hz, which is known to be a band spectrum, i.e. an absolutely continuous 
spectrum which is a union of closed intervals, which are called the bands. 



2.2.1 Floquet's theory 

We can first use the Floquet theory (or the Bloch theory, which is an alter- 
native name for the same theory). This is more detailed in the appendix. 
One can show that the spectrum of Hz is obtained by taking the closure of 

^k<^[o, 27V /T]cr{Hz,k) where 

is considered as an operator on L^(]R/TZ). So 

^(Hz) = Ufc,[o,t]^(^.,fc) • (2-7) 

We now write 

r = TZ and T* = — Z . (2.8) 

Hence we have to analyze for each k the operator Hz^k on L^(R/r). Later 
we will use the notation 

Hr = Hz,o. (2.9) 

A unitary equivalent presentation of this approach consists in analyzing 
Hz restricted to the subspace J)^ of the u G Lf^^{R) such that 

u{z + T) = e'''^ u{z) . (2.10) 

Here we did not see a fc-dependence in the differential operator but this is 
the choice of the space i)k (which is NOT in L^(R)), which gives the k- 
dependence. Condition fl2.10p is called a Floquet condition. 
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This means that we have written, using the language of the Hilbertian- 
integrals, the decomposition 

/■e 

L2(R) = / i)f^dk (2.11) 

J[0,2tt/T] 

and that we have for the operator the corresponding decomposition 

H,= r H^^kdk, (2.12) 

J[0,27r/T] 

with Hz^k unitary equivalent to H^^k- 



For each k G [0,27r/T[, H^^k has a discrete spectrum which can be de- 
scribed by an increasing sequence of eigenvalues {Xj{k))j^f^. The spectrum of 
Hz is then a union of bands Bj, each band being described by the range of Xj. 
At least when we have the additional symmetry even, one can determine 
for which value of k the ends of the band Bj are obtained. For j = 1, we 
know in addition from the diamagnetic inequality that the minimum of Ai is 
obtained for k = : 

inf Ai(A;) = Ai(0). (2.13) 

k 

2.2.2 Wannier's approach 

When the band is simple (and this will be the case for the lowest band in 
the regime e small), one can associate to Xj{k) a normalizecH eigenfunction 
ipj {z, k) with in addition an analyticity with respect to k together with the 
(27r/T)-periodicity in k. 

In this case (we now take j = 1), one can associate to ipi, which satisfies, 

ip,iz + T;k)=ip,iz,k), (2.14) 

and 

ipi{z;k + y) = ¥>i{z,k) , (2.15) 
a family of Wannier's functions {ipi)i(zr defined by 

tpo{z) = — I exp{ikz) Lpi{z, k) dk , VK^) = - i) , (2.16) 
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for £ e r . 

In addition, we can take ipo real. One can indeed construct ipi satisfying in 
addition the condition 

ipi{z, k) = ipi{z, -k) . (2.17) 
One obtains (after some normalization of ipo) that 
Proposition 2.1. 

(i) The family {ip^)f^Y gives an orthonormal basis of the spectral space at- 
tached to the first hand. 

(a) ipQ is an exponentially decreasing function. 

The second point can be proved using the analyticit}{§ with respect to k. 
This orthonormal basis corresponding to the first band plays the role of the 
basis Pj{z) exp — in the Lowest Landau Level approximation. 
Note that we recover ipi{z, k) by the formula 

ipi{z,k) = exp{—ikz) ^^exp(iH) ipi{z) . (2.18) 

Moreover, the operator A on i'^{T) whose matrix is given by 

Aw = {H,iJe,i^e') (2.19) 

is unitary equivalent to the restriction of to the spectral space attached 
to the first band. 

One can of course observe that A commutes with the translation on ^^(F), 
so it is a convolution operator by a sequence a E (F) (actually in the space 
of the rapidly decreasing sequences 5(F)), 

Aw = a{l - i') , (2.20) 

which is actually the Fourier series of t— > Ai(A;) 

Ai = a, (2.21) 

where 

Ai(^):=— / e-XY>{-ilk) \i{k) dk . (2.22) 



27r _ 
So we have 

{Au){l) = Y^a{t- i')u{i') , for M G ^^(F) 



''er 



^One can make a contour deformation in the integral defining ■00 in (|2.16p . 
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2.2.3 (iVT)-periodic problem 



There is another way to proceed at least heuristically. We keep w T-periodic 
but look at the (A^T)-periodic problem and we analyze this problem. The 
spectrum is discrete but the idea is that we will recover the band spectrum 
in the limit +00. If we compare with what we do in the Floquet 

theory, the analysis of the (A^T)-periodic problem consists in considering 
the direct sum of the problems with a Floquet condition corresponding to 

' NT ' ' NT 

Note that this decomposition into a direct sum works only for linear 
problems, so it will be interesting to explore this approach for the non linear 
problem. 

In this spirit, it can be useful to have an adapted orthonormal basis of 
the spectral space attached to the first N eigenvalues of the iVT-periodic 
problem (which can be identified with the vector space generated by the 
eigenf unctions corresponding to the N Floquet eigenvalues associated with 

' NT ' ' NT 

Our claim is that there exists an orthonormal basis, for the L^-norm on 
] — consisting of (A^T)-periodic functions and replacing the Wannier 

functions. 



We write 

1 ^ 

where is an eigenfunctioiill of the (A^T)-periodic problem, chosen in such 
a way that 

0f(z + T)=<-Vf(^), (2.24) 

with = exp{2i7i/N). 
We can then introduce 

= T/{NTZ), (2.25) 
and define, for £ G = F, the (A^T)- Wannier functions 

ij^{z)=^^{z-i) (2.26) 

^Note that except in the case j = 1, we do not claim that (f)j^ is the j-th eigenfunction 
but this is the first one corresponding to the condition (|2.24p . 
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This gives an orthonormal basis of the eigenspace attached to the first N 
eigenvalues of the (iVT)-periodic problem. These first N eigenvalues belong 
to the previously defined first band. 

Note that conversely, we can recover the eigenfunctions from the ■^j^ 
by a discrete Fourier transfrom. In particular we have 

j=o 

Except the fact that these "Wannier" functions are NOT exponentially 
decreasing at oo (they are by construction (A'"T) -periodic), one can then 
play with them in the same way (this corresponds to the replacement of the 
Fourier series by the finite dimensional one). We then meet the "discrete 
convolution" on i'^ij?'^) : 

Of course £^(r^) is nothing else than with its natural Hermitian struc- 
ture. 

We have presented different techniques to determine the bottom of the 
spectrum of H^, which all provide the same ground energy. We will now recall 
more quantitative results based on the so-called semi-classical analysis. 

3 Semi-classical analysis for the T-periodic 
case 

3.1 Preliminary discussion 

Till now, we have not strongly used that we are in a semi-classical regime: 
our semi-classical parameter here will not be the Planck constant h (which 
was already assumed to be equal to 1) but e. We will now use this additional 
assumption for extracting quantitative results from the previously presented 
qualitative theory. As already said, the physics literature is analyzing a very 
particular model, the Mathieu equation. We will rapidly sketch how one can 
do this in full generality. For the one dimensional case which is considered 
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here, one can probably refer to Harrell |Haj (who uses techniques of ordinary 
differential equations) or to the book of Eastham [Easj , but we will describe 
a proof which is more general in spirit, which is not limited to the one dimen- 
sional situation (see Simon [Si], Helffer-Sjostrand |HeSjl| , Outassourt |Ouj ) 
and is described in the books of Helffer |He] or Dimassi-Sjostrand |DiSj| . 
As we have shown in the previous section, the good description of the first 
band, can be either obtained by a good approximation of Ai(A;) and ^i{z, k) 
as e ^ or by first finding a good approximation of the Wannier function 
introduced in (12.161) . which is expected to be exponentially localized in 
one well, or of the (A^T)-periodic Wannier function introduced in (12.231) . 

The analysis is done usually in two steps. First we localize roughly Ai(fc), 
then we analyze very accurately the variation of Xiik) — Ai(0). 
The first one will be obtained by a harmonic approximation and the second 
one by the analysis of the tunneling effect. 



3.2 The harmonic approximation 



We will provide the explanation in a general case containing the model con- 
sidered by Snoek [Sn] as a particular case. We recall that we work under 
Assumption 11.11 The statements below are sometimes written vaguely and 
we refer to |DiSj| or |Hej for more precise mathematical statements. 
For the approximation of Ai,z(0) (actually for any Xi^z{k)) the rule is that 
we replace w(z) (having in mind (11.71) ) by its quadratic approximation at 0. 
The harmonic approximation consists in first looking at the operator 

1 cP w"(0) 2 , , 

on M. 

For the model in |Sn], w(z) = sin^(^), and we find 

This operator is a harmonic oscillator whose spectrum is explicitly known. 
The j-th eigenvalue is given by 



Af; = :^ ^yV^. (3.3) 
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The two main pieces of information we have to keep in mind are that the 
ground state energy is 



AtT = ^v/^, (3.4) 

and that the gap between the first eigenvalue and the second value is given 
by 

^har yhar _ ^kar ^ 1 . (3.5) 

The corresponding positive normalized ground state is then given by 



2 

^ljh'^'{z) = 7r-3w"(0)5 e"^ exp -w"(0)^|- . (3.6) 

It will also be important later to have the computation of the norm. So 
we get by immediate computation : 

ij^'^'^iz)^ dz = vr- W"(0)^ e-5 . (3.7) 

The mathematical result is that this value provides a good approximation of 
-^1,^(0) (and hence of the bottom of the spectrum of H^) with an error which 
is 0(1) as e ^ : 

AM(O) = AtT + 0(l). (3.8) 

By working a little more, one can actually obtain a complete expansion of 
^■^1,2(0) in powers of e and hence, of e\i^z{k)^ since they have the same expan- 
sion. For each j G N*, one has a similar expansion for eAj^2(0). This implies 
in particular an estimate of A2,2(0) — Xi,ziS^)^ called the longitudinal gap : 



5. := A2,.(0) - Ai,,(0) = ^ ^ + 0{l) . (3.9) 



From now on, we simply write Ai,^ or Ai instead of Ai,z(0) for the ground 
state energy of the periodic problem. 

Let us note that the ground state of the harmonic oscillator also provides a 
good approximation of the ground state of Hl^"^ . So we obtain, using (13. 7p 
that for 01, the L^-normalized ground state of if^^^, we have 

/•+¥ 

/ ^ (t)i{zfdz = 7r-5w"(0)3 + 0(1) . (3.10) 
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3.3 The tunneling effect 



We now briefly explain the results about the length of the first band, which 
is exponentially small as e — > 0. The results can take the following form 
(see the work of Outassourt |Uuj or the book by Dimassi-Sjostrand, Formula 
(6.26)) 

Ai(A;) - Ai(0) = 2(1 - cos(A;T))r + 0(exp -^-^) (3.11) 

with a > (arbitrarily close from below to 1) and, for some Cr 7^ 0, 

3 5* 

r ~ e~2 exp . (3.12) 

Moreover one can express the constants and S once w is given (se^ also 
|Hej in addition to the previous references). This r seems to be called in 
some physical literature the hopping amplitude. 

Here, we simply explain how one computes S which determines the exponen- 
tial decay of r as e ^ 0. In any dimension, 5* is interpreted as the minimal 
Agmon distance between two different minima of the potential w. In one 
dimension, with w satisfying Assumption (11.11) . this distance is simply the 
Agmon distance between two consecutive minima and is given by 

S:=V2 I Vw(;z) dz. (3.13) 



2 



In particular, when w(2;) = sin (^), we get 



S:=V2 \sm{—)\dz=^ . (3.14) 

I T 1 IT 

This is to compare to (14) in |SnStlj . which is not an exact formula (as 
wrongly claimed) but only an asymptotically correct formula. It can be 
found, for this Mathieu operator, in |AS] . 

Let us give the formula for the constant Cr- It can be found in [Ha], see also 
[Uu] . Formula (4.14) and ^ p. 58-59. We have : 

3 _i 

Cr = 2^'K aexpA^, (3.15) 

with (assuming w even) 

Ar = \im\ / / dz+ ^ In?? I . (3.16) 




^The computation is a little simpler in the case when w is even. 
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We just sketch the mathematical proof. Filhng out all the wells suitably 
except one (say 0), we get a new potential w^""^ > w which coincides with w 
in an interval containing and excluding small neighborhoods of all the other 
minima. We consider, for e small enough, the ground state of this modified 
problem and (multiplying by a cut-off function) we get a function ip^^ (and 
an eigenvalue A'^'^'') which is a very good approximation of_0o- 
Now the hopping amplitude in the abstract theory is giveiij exactly by 

- r = a(T) = {H^^o , ^i) = {{H, - /i)V'o , ^i) , (3.17) 

the last equality being satisfied, due to the orthogonality of ipo and ipi, for any 
fi. When replacing ipo by its approximation, one has to be careful, because 
tpQ^^ and := ipo^^i- — T) are no more orthogonal. So this leads to take 
/i = A^^^, and one can prove that 

r~-((i7,-Ar)^o"'', ^D- (3.18) 
An easy way to see that r is exponentially small is to observe that 

{{H, - , V^D = ((w(^) - w^^'^W^ , V'D , (3.19) 

and to use the information on the asymptotic decay of ip'^^ . The WKB- 
approximation of ip'o^ is, in a neighborhood of 0, 

= e-3 6(2,e)exp — / ^w{s)ds, for ^ > , (3.20) 

with 

6(z,e)~5^6,(z)e^ (3.21) 

i>o 

and 




(3.22) 



It should then be completed by symmetry to get an even WKB solution on 
]-T,+T[. 

Note that we have 




which implies that 6o tends to +oo as 2; ^ r_. 



^For the Mathieu potential, this is consistent with Formula (13) in [SnStl] , 
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An integration by parts together with a WKB approximation leads to the 
asymptotic estimate of r announced in fl3.12p . More precisely, we get that 

the prefactor is immediately related to the constant 6o(|-)^'y/^(|-) and 
this leads to (I3.15p . Note that more generally we have 

boizMT - z),/Vizj = Cst, (3.23) 

which again shows the blowing up of 60 at T. 



Finally, we emphasize that i/jq is a good approximation of i/jq only in 
intervals ] — T + ri,T — ri[ for some rj > 0. 

One can also see that a{kT) is of the order of |a(T)|l'^l (for k > 2) 

a{kT) = d{T^) , (3.24) 

so it is legitimate in order to compute the width of the first band to forget 
allthea(£) ioi i eT,£ ^ 0,±T. 

Thus, in the k variable, the spectrum (corresponding to the first band) is up 
to a very small error, of the order of the square of a(T), given by the operator 
of multiplication in L^(]R/r) by the function a(0) + 2a{T) cos(/cT). 



Remark 3.1. 

What is written above corresponds to the use of Wannier functions on R. One 
can write a close theory using the (NT) -periodic Wannier functions without 
modifying the main terms of the asymptotics. In particular, ip^^^ is also a 
good approximation of ipo in intervals ] — T + 77, T — ?][ for some t] > 0. 
The interest of the Wannier functions on M is that they allow to recover the 
information for all Floquet eigenvalues (see the discussion in Section\7J\). 



4 Justification of the reduction to the longi- 
tudinal energy 

4.1 Main result 

In this section, we address the reduction to the energy defined in f ll.35p 
and prove the following theorem (recall that is defined in f ll.Sip ): 
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Theorem 4.1. // 

(Ana) m^{e,g){uj^-n)~^ «1 (4.1) 

and 

(Anb) g{2tu±-n)m'^{€,g)(ij±-n)-^ «1, (4.2) 

wc have 

inf Sr'^m = u^± + ?)(1 + 0(1)). (4.3) 

11*11=1 

Both Theorem 11.31 and Theorem 11.41 are a consequence of Theorem 14.11 as 
soon as we have the appropriate rough estimates on aheady presented 
in the introduction. This is what we explain first in Subsection 14.21 before 
proving the theorem in Subsection 14. 3[ 



4.2 Proof of Theorem 11.31 and Theorem 11.41 

4.2.1 Weak Interaction case 

In the Weak Interaction case, we recall from f ll.57p . that, when (11.550 is 
satisfied, then 

^ 1/e. (4.4) 

Therefore, when f ll.54p and fll.55p are satisfied, then (14. ip and (14.20 auto- 
matically hold with the observation that 

g{2uj^-n){uj^-nY-^m^{e,g) < C g{2uj^ - n)e^{{uj^ - n)e)-^^ « 1, 

and Theorem 11.31 follows from Theorem 14.11 

4.2.2 Thomas-Fermi case 

In the Thomas- Fermi case, we will prove in (15.180 that, when (11.590 and 
(ll.60p are satisfied, then 

ml ^ {gujjef/' . (4.5) 

Let us verify that, if ffL59D . ffLHCiD and ffLHTl) are satisfied, then (O) and 
(14.20 hold. This will prove Theorem 11.41 

We get (14. ip in the following way. First we have : 

{u^ - a)-'m%,g) < C{u;^ - Qy'ujgh-"^ . 
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Hence (14. II) is a consequence of 



guj±_ « e(a;_L -^)^ , (4.6) 

which follows from fll.6ip since (11.591) and (11.611) imply that {uo^ —VL)t » 1. 
The check of (14. 2 p is then immediate from fll.6ip and fl4.5p . 



4.3 Proof of Theorem 14.11 

Because of the upper bound (ll.53p . Theorem 14.11 is a consequence of the 
following proposition, recalling that 5± = uj± — VL . 

Proposition 4.2. 

There exists a constant C > such that, for all e g]0, 1], for all uj^, VL s.t. 
(5_L > 1 and for all g > , 

M Q'-'''{^)=uj^+m^{e,g){l-CrA{e,g)) , (4.7) 
11*11=1 

with 

1 

< r^(e,?) < g'/'6-J (^^^) ' rn^{e,g)i + mA{e,g)6^-' . (4.8) 
Proof of the proposition 

For simplicity, we make the proof for = 0. Note also that 

l-CrAie,g)>0 

by the lower bound. So we have only to prove (14. 8 p under the additional 
condition that the right hand side of (14. 8 p is less than some fixed ao . In any 
case, the estimate is only interesting in this case ! 

The proof does not depend on and for Q not zero, we will make a remark 
at the end on how to adapt it, using the diamagnetic inequality. 

The proof is inspired by |AB] where a reduction is made from a 3D to a 2D 
setting for a fast rotation. We project a minimizer onto tp± ® L'^{M/NT7j), 
and call ip±{x, y) ^{z) its projection: 

'^{x,y,z) = il)±{x,y)^{z) + w{x,y,z) with I %jj±{x,y)w{x,y, z) dxdy = . 

(4.9) 
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The orthogonality condition imphes in particular 



. iVT 

1= / ' \^(^z)\'^dz + / \w{x,y,z)\'^dxdydz (4.10) 



and we have the lower bound 



NT 
2 



NT JM2x]-i^,i^ 



with 



£'q{w{-, z)) dz > {6± + uj±) I \w{x,y, z)]"^ dxdydz , 

(4.11) 

S'sm = Q|V.,,^(a:,y)r + ^ix' + y') mx,y)\'^ dxdy . 

We compute the energy of \I' and use the orthogonahty condition and the 
equation satisfied by ip± to find that all the cross terms disappear so that 

NT 

Q^.^'-(vl/) = [ ' \^{z)\' dz + iO 

I NT 

NT 

+ ^a' iu){x,y,-)) dxdy + £'^{w{-, ■, z)) dz 

+ g \'^{x,y,z)\'^dxdydz, (4.12) 



where 



. NT 

^a'(0)= / 1 (l\<P'iz)\' + W,{zM'] dz. 



2 



NT 
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From flCTj) . KTTh and flCTj) . we find 



NT 



g^'P-(vl/)>^^ + ^^ /■ ' S's{w{;;z))dz+ I Sf{w{x,y,-))dxdy. 

d±+UJj_J_NT J^2 

(4.13) 

We use (14.131) together with the upper bound (11.531) and (14. lip to derive that 

\w{x,y,z)\^dxdydz<'^^^^^. (4.14) 



Note that the righthand side in (I4.14p is very small according to Conditions 

gl]) and m. 

Note that KIM implies 



\m'dz>l-^^^^^. (4.15) 



NT 
2 
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Then, we get also, 

I V.,.^(x, y, z)\' dxdydz < 2^ ^^^^^^ 
Ir2,,]_nt nt^ \dzw{x,y, z)]"^ dxdydz < 2 m^{e,g) . 

The proof of the Sobolev embedding of if^(R'^) in L^(M^) gives (see for ex- 
ample |Brej . p. 164, line -1) for a general function v in H^{E?) 

\\vh<^^M\T\\^yv\\y'\\^M\y^■ (4.17) 
Here || ■ ||p denotes the norm in Lp(M'^). 

In our case, we are working in H^{E?^y x (M^/A^TZ)). A partition of unity 
in the z variable allows us to extend this estimate also this case, and we get, 
for another universal constant C, 

Iklle < CN\\dM\T\\dyw\\T {\\dM\l + Ikll')'^' , (4-18) 

where this time || ■ ||p denotes the norm in Lp(M^ ,^x] — ^[). 
So we obtain : 



\\wh<Cm'l{e,g)-^[-^j-^j . (4.19) 

(C, C are A^-dependent constants possibly changing from line to line.) 
Since by Holder's Inequality, 

II II ^ II l|l/4|| ||3/4 



we deduce that 



\w\U < C mAie,g)Ur'^ (^^f^) ' " ^^"^O) 



We expand 



i^r = i^±iier+2i^±nei>p+4(3ft(^xe^)+^kp)'+4i^^nep3fj(^^e^) , 

Since f l4.12p implies that 

S''{^)>u;^ + S^{0-'^9 I \iP^{x,y)W{z)\^\w{x,y,z)\dxdydz, 
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in order to get the lower bound, we just need to prove that the last term is 

a perturbation to 

We can do the following estimates 

9 J \i^±ix,y)f\^{z)f\w{x,y,z)\dxdydz 

< Cogu±^J \^p±{x,y)\*dxdy)^ (/ \^{z)\^dz)^ \\w\\4 

< c,g'/\S^{Om\w\U 

< C2g'/'S^~'^ (^)' m^{e,g)Hs^{Or/' 

<Csg'/'Sr'^ m^{e,g)-. {l + C m^{e,g)Sr') £^(0- 

Here to get the last line, we have used the lower bound 

and Km . 
This leads to 

S^{^)>uj^+S^{0 (^l-Cg'/'Sr's ^^^^y m^{e,g)l-Cm^{e,g)S, 
and then to K7} . 



Remark 4.3. 

In the case with rotation Q, the proof is the same if we replace S'^ by S'^ q 
defined by 

^knW = j^^ (^|V.,,^ - inr^i^? + - fi')rVl') dxdy . (4.21) 

We also use the diamagnetic inequality 

J \\/\w\{x,y)\^dxdy < j \ {Vw - iVtr^w) {x,y)\^ dxdy (4.22) 

which provides the Sobolev injections. 
Remark 4.4. 

Here, we have not proved that the minimizer of S behaves almost like the 
ground state in x,y times a function of ^ which minimizes £a- We are just 
able (see fl4.14p ) to prove that the minimizer is close to its projection (in 
some or norm). When N = 1, this can be improved under the stronger 
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condition (11 .491) . We first observe (note that fl4.13p is still true with the 
addition of S'^{C,) on the right hand side) that 

S'^{0<mA{e,g). (4.23) 

Using (14.151) . assuming < 1, this leads to 

£m<mAie,g)a-'^^^^)-m' (4-24) 

We will show in Subsection \5.^ (see (15.141) ) how to proceed in order to show 

that ^ is close to the ground state (f)i{z) of H^^^ . 

This can allow to improve the information given in Theorem \1.2[ 



5 The ID periodic model : estimates for 

The aim of this section is to analyze m^- We note that rough estimates 
were aheady given for the weak interaction case which were enough for the 
justification of the model but the corresponding rough estimates needed for 
the Thomas-Fermi justification will be obtained in this section. We will then 
look at accurate estimates for m^, which will be established under stronger 
hypotheses. We will end the section by the discussion of the case > 1, 
which finally leads to the introduction of the DNLS model for the Weak 
Interaction case. 



5.1 Universal estimates 

We consider the one dimensional situation and a T- periodic potential W, 
which could be for example W{z) = (sin7r2;)^/e^. We consider the problem 
of minimizing on L^(]R/r]R) the functional 

T T T 

V^^^W = ^ W{z)\'dz+ r W{zmz)\'dz + g r \i^{z)\'dz, 

2J-T J_T J_T 

(5.1) 

over \\iP\\l2 = 1. 

We are interested in the control of the minimum of the functional and will 
simply prove 
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Lemma 5.1. 

Ifg >0, then 



m^^j:= inf G (i;) ^ X, + g [ " dz + 0(9) , (5.2) 



where (Ai,0i) is the spectral pair of — | ^ + W{z) corresponding to the 
ground state energy (with ||0i|P — i)- 



Proof : 

It is clear that 



Ai <m(^) < Ai+^ / 101(^)1^0?^, (5.3) 



so the question is now to improve the lower bound. 

One could of course think of applying bifurcation theory but this gives only 
a local result and we need in any case a global estimate for showing that the 
global minimizer of Q is closed to (pi as 'g is small. 
Let (pmin be a minimizer of Q, then we know that 

T T T 

dz+ / ' W{z)\<Pmin{z)\'dz<Xi + g I' \M^)\Uz. (5.4) 

It It 

2 ''-2 ''-2 

So 0mj„ plays the role of a quasimode (or approximate eigenfunction) for 

A rather standard theorem in perturbation theory (we can write <^rain — 
CK^i + gives first 



1 /" 2" 

2 / ^ \^rmnf^ 



1- = ll^^ll^ <g. 



l\\<\>i{z)\Uz 

2 

A2 — Ai 

and then the existence of a complex number c of modulus 1 such that 



T 



r\\uz)\'dz 

||0m^n " C^iHi^ < 2^^- ^ . (5.5) 

Here A2 denotes the second eigenvalue of or Hamiltonian. 
Of course, the estimate is only interesting if 

f\_\<^Az)\^dz 
2? ~~\^ _ < 1 ■ (5.6) 
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We can now write 



fn{g) > \i+g |0mm(^)|^ dz 



> Al \(j)i{z)\^ dz - 4:g\\(pmin - C0l||2||0l||6 

2 

For the last estimate, we develop l^mml^ in the following way 



(5.7) 



'mtn 



1^ = |C01 + (j)min - C0l|^ 



(5i 



> |0i|^ + 2|0i|2|0^i„ - C0ip - 4|0i|3|0^i„ - C0i| . 

From this inequality we get 

> |0i|^-4|0in0^„,-c0i|. (5.9) 



■'mm I 



It just remains to control ||0i||6 uniformly with respect to ^, which can be 
deduced of the uniform control of the norm of 0i in L^. ■ 

One can actually be more precise on what we have claimed in Lemma ISTTl 

Lemma 5.2. 

//^ > 0, then 

m{g) > \i+gm\\ - 2igi\\ct>i\\l\\ct>i\\l{\2 - Ai)"^ . (5.10) 
This estimate is interesting if 

?<^ll0i|l4ll0ill6'(A2-Ai). (5.11) 

Remark 5.3. 

Everything being universal, one can of course replace T by NT in the de- 
scription. 

5.2 Semi-classical results in the Weak Interaction case : 

N = 1 

We first recall that using ( IB.lOp we have, under Condition (11.551) . the rough 
control 

T 

1 , . . . , ^ n . , c 



Ce 



<Kz<mA{e,g)<\z + g |0i(z)|^rf2;< -, (5.12) 
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which leads to fll.57p for = 1 and was sufficient for the justification of the 
longitudinal model A. 

Let us now show that under stronger assumptions one can have a more 
accurate asymptotics including the main contribution of the non-linear in- 
teraction. 

Proposition 5.4. 

Under Assumption (11.491) . uia admits the following asymptotics : 

m^(e, g) = Af ^(e) + 7r-^w"(0)^?e-5 + cq + 0(e) + 0(?ie-^) . (5.13) 

Proof : 

Indeed, Ai and Ai — A2 are of order i, and by (I3.1UI) and (15. 5p . we get 

W'Pmin - < Cge^ . (5.14) 

Using the harmonic approximation, the term ||0i||6 is of order e~s and the 

3 1 

remainder appearing in (15.101) is of order ^2 6~4. Altogether we get for the 
energy 

mAie,g) = X,,, + g I \Mz)\Uz + 0{gh--^) . (5.15) 

Using (I3.10p . we obtain (15.130 . This asymptotics becomes interesting in the 
semi-classical regime if (I1.49P holds. ■ 

Remark 5.5. 

Exponentially small effects will be discussed in Section^ 

5.3 Semi-classical analysis in a Thomas- Fermi regime : 
case = 1. 

5.3.1 Main results 

In this subsection, we ffist give the rough estimate leading to (ll.62p for A^ = 1. 
Recall that g = 9^±, but 'g and e are taken as independent parameters. 

Proposition 5.6. 

If for some c> 0, 

ge^<c, (5.16) 
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and if 

ge^»l, (5.17) 
then there exist C and eo such that 

^gie-i < mA{e,g) < Cgie'i , Ve e]0,eo] . (5.18) 

We will also get the following accurate estimate : 
Proposition 5.7. 

ge'KKl, (5.19) 

and (15.171) are satisfied , then 

mAie,g) = 2-1335" V(0)%3e-i (^1 + 0{g~h~^^)^ . (5.20) 

The new asusmption is (15.191) . which is stronger than (15.161) . 
5.3.2 The harmonic functional on R 

2 

Let us start with the case of a harmonic potential We{z) = 7^ on M, with 
7 > 0, and consider the problem of minimizing 

z n T 



g^-.^ (u) = - I u\tY dt + ^ I ^ t'uity dt + g I ^ U(tf dt (5.21) 



T It 



over the m's in the form domain of g^^'^ such that = 1. 
We denote by m^^''^ the infimum of the functional. Actually there are two 
approximating " harmonic " functionals of interest corresponding to T finite 
and to T = +00. An interesting point is that, for T large enough, the 
minimizers of these two functionals are the same as we will see below. But 
let us start with the case T = +cxo. 

Lemma 5.8. 

If dSHD holds, then 



Hr,+oo 



'e,g) = 2-hh~^-f^gh-^^ (l + 0{g-h-^s)^ . (5.22) 
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The proof is very standard (see for example |BBHj . |Af] or |CorR-DY] 
which treat the (2D)-case). The analysis is done through a dilation. We 
look for an L^-normalized test function 6 in the form 



[z) = p-^v{pz) , (5.23) 



with p and v to be determined. 
The 1 — D energy of becomes 



2 

with 



^ p" I v\tfdt + p-^e-^ / ev{tfdt + gp / v{tYdt, (5.24) 



e7 = e/\/^7 




This leads to choose p = p-y such that 



hence 



p^ = e^^^3, (5.25) 
and the energy of this model becomes 



?3e-3 [qTF{v) + -(e|^)"3 j^v'{tfdtj (5.26) 

with 

qTF{v):= / t^v{tfdt+ I v{tf dt . (5.27) 
Jr Jr 

2 2 

This is asymptotically of the order of ^3e~3 and Condition fl5.17p is just 
the condition that the kinetic term is negligeable in the computation of the 
energy. 

Let us recall the details of this asymptotics for completeness. We have first 
to minimize over L^-normalized v the approximating functional qte- The 
minimizer Vmin(t) of qte is determined by the equations 



rv{t) + 2v{tY = Xv{t) , / v{tydt = 1 . (5.28) 

We get 



1 

^2 



Vrmnit)=2-HX-t')l, (5.29) 
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with 




(5.30) 



and for x G M, 

{x)+ = max(x, 0) . 
The corresponding TF-energy is 

CtF := / {t^VrrUnitf + Vmin{t)^) dt = . (5.31) 

Unfortunately, this minimizer is not in H^{M.) and can not be used directly 
for our initial rescaled functional 

Qtf{^) = <1tf{v) + (r v'{tf dt , 
JR 

with 

a = g -it^K 
Here we recall that (15.171) implies 

< a « 1. 

So we need to regularize this minimizer to have an upperbound for the energy 
of our "harmonic" functional which is good as cr 0. 



This can be done in the following way (see for example [Af] and references 
therein). 

We introduce ^ 

lis) = i -1 -r / i (5-32) 
{ sa 6 , if s < crs 

Let us consider the function 

We get that Vo- belongs to and satisfies 

j \v^{t)\'' dt = I - r{a) , 

with 

r{a) = 0{J). 
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More precisely the (positive) remainder r(cr) is 

r{a) = / ^{-\Vniin{t)\'^ + o-~^v.m.in{t)\^)dt. 

Let us now consider test function 
Then we have 

1 4 

= (l + -r((T) + C((T3)Xv. 

So we get 

a 

So we have the upper-bound in the statement (15.221) (actually with a 
better remainder term) of the lemma. The lower bound in (15.221) is immediate 
because the kinetic term is positive. 



5.3.3 The harmonic functional on ] — 

We consider now the case of the interval and have the following Lemma : 
Lemma 5.9. 

Under Assumption (15.170 . there exists C > such that 

m'r^{e,g)>^gh''^. (5.33) 

The proof is based on the same method as in the previous subsubsection. 
It is easy to see that the minimizers coincide if 

^>\K (5.34) 

that is 

1 ^ /3\ ^ 

T>g-^e^\^-j . (5.35) 

If (I5.35P is not satisfied, we can still have a lower bound for the infimum 
of the functional. The renormalized functional reads 

pT pT pT 

qren,T^^^ ;= p2 j' ^ ^' ^t + p~\:^^ j^' t^v{tfdt + gp v{tfdt, (5.36) 
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pT ' ■ - 7 ■ • / pT 

2 2 2 



which satisfies 
Using the Holder inequahty, we obtain, if ||t>||2 = 1, 
and using our assumption, we obtain 



g---T^^)>lA-^(^p)>l^le-i, (5.37) 

if ||t;||2 = 1. 

We then immediatly obtain Lemma I5.9[ 

5.3.4 Relevance of the "harmonic functional" for rough bounds 

First we prove Proposition 15.61 We can proceed by direct comparison. Ob- 
serving that we can find a > such that 

^{z)<az\ V2;e [-|,+|], 

and 

p^T > 2A5 . 

Here, we use (15.161) and 

1 , _ 1_ 1 _ 1. 

Pa = Coa3(e »g ») > Coa^c 3. 

We can then use the asymptotic estimate (15.221) with 7 = a to get the upper 
bound in (I5.18p . 

Using now Assumption (11.11) . we can also find a such that 
wiz)>az\ V2;G [-|,+|], 
This leads, using our analysis of q^^ in the harmonic case to the lower bound 

in (Km> . 
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5.3.5 Relevance of the "harmonic functional' for the asymptotic 
behavior 



In order to have a better localized minimizer, we should assume that p 
+00 and this corresponds to replacing Assumption (15.161) by the stronger 
Assumption fl5.19p . 

Moreover, we have to verify that under this assumption the "harmonic 
approximation" is valid for this energy computation. For this, we should 
analyze the localization of the minimizer. Assuming that such a localized 
minimizer exists (minimize the functional v ^ j{z'^v{zY + v{zY)dz)^ we 
can also get an upperbound of by using a harmonic approximation and 
a lower bound of the same order. 

For the lower bound, we have just to analyze (forgetting the positive 
kinetic term) the infimum of the functional 



T 

2 / W[Z: ,2 



+ I dz . 



2 ^ 

As in the other case, a minimizer (over the L^-normalized (/''s), should satisfy, 
for some /i > 0, the Euler-Lagrange equation 

+ 2gct>{zf = fi<f){z) , 

where /i will be determined by the normalization over ] — ^, ^[. 
We find 

A-^y. (5.38) 



V29 

with 



e2 



— / {fi ^)+dz = 1. (5.39) 



1 f w(z). 

But we know from the upperbound that /i is less than two times the energy 
which is asymptotically lower than m^^leg). In particular, if /ie^ is small, it 
is easy to estimate /i using the harmonic approximation of w at its minimum. 
It remains to verify the behavior of /ie^. We find 

Not surprisingly, this shows that /xe^ is small as p — > +00. So finally, we have 
obtained Proposition 15. 7[ 
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5.4 The case iV > 1 



We would like to extend our rough or accurate estimates for to the 
analogous estimates for > 1, keeping the same kind of assumptions. 



5.4.1 Universal control 

We now consider the functional over ] — Using the minimizer ob- 

tained for = 1 and extending it by periodicity, we get after renormalization, 
the general upper-bound 

m^(e,?) <m^(e,^). (5.40) 

From this comparison, we obtain immediately the rough upper bounds in 
the WI case and in the TF case. 



5.4.2 Rough lower bounds 

In the Wl case, we always have, observing that Ai^^ is the ground state energy 
for any N eW, 

Xl<m^{e,g). (5.41) 
Hence we obtain in full generality 
Proposition 5.10. 

Under Condition ([i.5^[j, then, for any N > 1, we have 



m^{e,g)^-^ (5.42) 

In the TF case, it remains to prove the lower bound which will be a 
consequence of the following inequality : 

m^{e,g)>^ghl. (5.43) 



We indeed observe that if mat is a normalized minimizer, then there exists 
i,(j + 2)f 



one interval Ij :=]jf , (j + 2)f [ (j G {-A^, . . . , - 2}), such that 



\uN\'^dz > 
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We can then write, forgetting the kinetic term and translating Ij to ] — -j, +|-[, 

m;^(e, 'g) > e"^ w(2;) \un\'^ dz + 'g Imat]^ 

> inf(||MAr||2, IImtvII^) inf||„||=i {W,\u\'^ + g\u\'^) dz . 

2 

Then we can combine the lower bound obtained for = 1 and the inequality 
w(z) > az'^ to get f l5.43p . So we get finally that has the right order in 
the TF case. 

Proposition 5.11. 

Under Assumptions (15.161) and (15.171) . we have, for any N > 1, 

m^{e,g) ^ g^e^ . (5.44) 

This extends to general our former Proposition 15.61 
5.4.3 Asymptotics 

We would like to give conditions under which the universal upperbound (I5.40p 
becomes actually asymptotically or exactly a lower bound. 

Proposition 5.12. 

Under either Assumption (I1.49P or Assumptions (I5.17P and (I5.19p . 

"^A (e, g) ~ mA{e, -^) . (5.45) 

Proof : 

The upperbound was already obtained in (15.400 . The proof of the lower 
bound is different in the two considered cases. 



WI case. We will see later (in (17. 6p ) by a rough analysis of the tunneling 
effect and the property that the infimum of the function 

N-l 

3 (Co, C2, . . . , Cn-i) ^ |Cj|^ 

i=o 
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over ^ . |cjp = 1 is attained when all the |cj|'s are equal 



|c,| = , for j = 0, . . . , AT - 1 , (5.46) 



that, under Assumption (11 .49^ . there exist C > 0, cq > and a > such 
that 

m^((?,e) >mA(^,e)-C(?+l)exp--, Ve G (0, eo] . (5.47) 

IS e 



TF case. In this case we can for the lower bound forget the kinetic term 
and come back to the analysis of Subsubsection 15.3.51 with T replaced by 
NT. Under Assumption (15.191) . we have seen in (I5.38P that the minimizer 
mtv is localized in the neighborhood of each minimum and T-periodic. 
We can then write 

NT T 

J_It(^\un\'^ + g\uN\^) dz = N J^Ti^luNl"^ + g\uN\^) dz 

2 T ^ 

= Ht Iv^MjvP + ll ViVwivl") dz 
>inf|H|=i/_\ (^^|^|2 + ||^|4j dz. 

But under Assumptions (15.171) and (I5.19P , the last term in the inequality has 
same asymptotics as mA(e, ■§) and we are done. 



Remark 5.13. 

This proposition leaves open the question of the equality in (15.451) . 



6 Study of Case (B) : Justification of the 
transverse reduced model 

6.1 Main result 

We have defined S^q by (ll.40l) -( |1.41l) and ttl^q, the infimum of the energy 
by (11.521) . In case B, the proof of the reduction does not depend on whether 
A^ = 1 or A^ > 1. The only difference is when looking at the rough or accurate 
estimates of the reduced model. Note that only rough estimates are used in 
the part concerning the justification of the model. 

The reduction is very similar to case A, and we will prove 
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Theorem 6.1. 

[RBa] eml^«l, (6.1) 

and 

(RBb) gm%e"^«l, (6.2) 

then, as e tends to 0, 

„ inf gr^'l^) = Ai,. + ml^il + 0(1)) . (6.3) 
11*11=1 

Then Theorems 11.51 and 11.61 follow from this result and appropriate esti- 
mates on m^Q, as we will prove in section [U?^ while the proof of Theorem 16. II 
is made in Section 16.31 

6.2 Proof of Theorems 11.51 and 11.61 

The issue is to determine the magnitude of the infimum of the energy of the 
transverse problem m^^. 

6.2.1 Reduction to the case = 1 

As in Case A it is immediate to see that 

< "^B,^7(^,w±) • (6.4) 

If indeed ipmin,N was the T-periodic minimizer for (11.381) with = j^, we 
get ([631) by using p.26p, i\2.'27\i and taking i/jj^^^ = -^i)min,N- 



So it remains for the needed upper-bound to analyze the case N = 1. 
This depends on the magnitude of g and leads us to consider two cases. 



6.2.2 The Weak Interaction regime : case = 1 

Proposition 6.2. 

// M.6b]) holds, then 

mB,n(g, t^±) < Cu± . (6.5) 
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Indeed, f ll.66p implies that g is bounded and the test function (which 
is independent of VL) imphes the proposition. 

Therefore, if (11.661) and (I1.67P are satisfied, then Theorem 16. II holds and 
implies Theorem 11.51 

6.2.3 The Thomas Fermi regime : case = 1 

We start with the case when f2 = 0. When g is not bounded, we can meet a 
Thomas-Fermi situation. 

Proposition 6.3. 

Ifg — s> +00, the function mB(g,ujj_) satisfies 

'mB{g,uJ±) ~ CTF^^±Vd, (6.6) 

with 

ctf = = 3-^2^-^ . (6.7) 

Therefore, if (Il.TOp . (11.710 . (I1.72p are satisfied, then Theorem 16. II implies 
Theorem 11.61 

Proof. 



A rescaling in y y g/uj± yields a new energy 



[ [ ^\Vu\^ + ^/gr^\u\'^ + 2^/g\u\'^] dxdy . 



" 2 

which is of the type Thomas Fermi (that is kinetic energy can be neglected) 
if 

^ « Vd- (6.8) 
V9 

This leads then simply to the TF reduced functional 

M I— >• (cuj, a/^) / ( o^^l'^^P + l""!^ ] dxdy , 

whose infimum over the unit ball in L^(R^) is of order ctf{^±\/9)i with 
ctf > defined by : 

ctf= inf / (^r^\u{x,y)\'^ + \u{x,y)\A dxdy. (6.9) 
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The minimizer exists and is explicitly known as 

1 1 3 1 

Umin{x, 2/) = -(A - r^)4 with A = 2271^2 . 

This leads to (16.71) . 

In addition, by a computation similar to the one in Subsection 15.31 we 
obtain more precisely 

Lemma 6.4. 

There exists c such that, as 'g tends to +00, 

— = ctfV^ +^ In (? + 0(—=), (6.10) 

with ctf defined in (16.91) . 
Remark 6.5. 

Note that we have the universal lower bound 

mB(g,UJ±) > CTF^^lVd- (6.11) 

This lower bound becomes better than the universal lower bound by ujj_ as 
soon as 

CTFV^>1. (6.12) 

Remark 6.6. 

In the semi-classical regime, conditions (BTFa) and (BTFc) in Theorem \l.(A 
(take their product) imply that this two-dimensional energy is much smaller 
than 1/e, that is 

uj^g^e-^'^ «e'^ . (6.13) 

We now look at the case when ^2 > 0. The previous proof, using that the 
minimizer of the TF reduced functional in (16.91) is radial, yields 

Proposition 6.7. 

There exists C such that, asg—>- +00, 

mB,n{g,u}±) < mB{g,uj±) + C In gg'^ . (6.14) 

This will be improved in (I6.48p by a direct study of the minimizer of SB,n- 
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Remark 6.8. 

For a lower bound, we can use the TF reduced functional 
whose minimum is explicit : 



inf In{u) = u;±v^eTF\/ ^(1 - ^^V^±^) ■ 
ll"ll=i V ^ 

Thus we get that, if there exists (3 E [0, 1[ such that 

0<fi/cJx</3, (6.15) 

then, asg—^ +oo. 



mB,n{g,^^A.) ^ g . (6.16) 
The uniformity of the approximation depends on (3. 

In fact, if one wants a more precise expansion of the energy, one can 
use the ground state p of 1^ to split the energy SB,n{u). Indeed the Euler 
Lagrange equation for p multiplied by (1 — Iwp) for any function u yields the 
identity (see 

Sb,^ = In{p) + j P^\Vv - zfi X rt;p + ~gp\l - It;^)^ 

where v = u/ p. Thus, Iq always provides a lower bound with an inverted 
parabola profile as soon as we are in a TF situation. The second part of the 
energy has the vortex contribution which is of lower order when fl/LO± << 
1. More precisely, the first vortex is observed for a velocity Q of order 
uj±lng/y/g. When Q increases and becomes at most like (3uj±_ with (3 < 1, 
the two parts of the energy I{p) and the rest become of similar magnitude. 
In the limit, Q — >■ u±, there are a lot of vortices and the description can be 
made with the lowest Landau levels sets of states. The leading order term of 
the energy is the first eigenvalue of —(V — iQ x r)^ which is equal to Q. 



6.3 Proof of Theorem 16.11 

We recall that we have the universal upperbound fll.64p . The lower bound 
follows from the following proposition and the fact that there exists c > 
such that 

as e tends to 0. 
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Proposition 6.9. 

There exists a universal constant C > such that 

„ inf Qr^'m = Ai,. + <n (1 -Cr^)- (6.17) 
11*11=1 

with 

0<r^< <,ni€r'+9k€)-k<,n)kl + jfy^- (6.18) 

Before giving tlie detailed proof, let us shortly sketch the case = 1. 
The proof is indeed essentially the same as for Case {A). One has just to 
exchange the role of (A) and (B). ruA should be replaced by m^, by Ai^^ 
and 5z by 5^. Note that in the two models {A) and {B) the ratio between 
the ground state energy and the splitting is of order 1. We have indeed 

(5_L = uj±_ and, in the semi-classical regime, — ^ ~ 1 . 

Oz 



Proof : 

We start from a minimizer \l/ and first write 

$=niv^+w (6.19) 

where Hat is the projection relative to the first eigenf unctions of intro- 
duced in (II .331) . We have 

n^u; = 0, (6.20) 

and 

+ IIHtv^II^ = 1. (6.21) 

We have the lower bound 

£'a{w) dxdy > z / \w{x,y,z)\^dxdydz , (6.22) 

with 

NT , , 

E'M ■■= + dz . (6.23) 

We now rewrite the energy in the form 

NT 

QP--^^{m)= ' £'^^{m)dz^ £'a{Un^) dxdy+ I S'^{w)dxdy+lNm 



NT 
2 



(6.24) 
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with 

Ij^{^iJ) = g I l^Hl^dxdydz , (6.25) 

and 



I {\\V,^yij-inr^ij^ + Uu^^-n^y\'ipA dxdy, (6.26) 

■In., \^ ^ J 



with r_L = (— x). 

We note that /at > and that 

^^,nW> ^±11^11'- (6.27) 
We first get the control of Having in mind f ll.64p . we obtain 



and this implies 



M 1,9 > TM9 (6-28) 



N 

kir<^- (6.29) 



The right hand side in 06.291] is small according to (16.11) . Note also that we 
have immediately from (I6.2ip . 

N 

||n;v^|P>l-^. (6.30) 

We now have to control the derivatives of w. For the transverse control, we 
start from 

>^i,z + n > \z + t; {V^^yW - iQrj_w\'^dxdy , (6.31) 

2 J^2_^y^]_NTN 



'^x.y^i 2 ' 2 



which leads to 

\\\'^x,y'W — i^rxwlW"^ < 2m^ Q . (6.32) 
For the longitudinal control, we write, for any a G [0, 1] 

Kz + m^,, > Ai,,||n;v^||2 + |||9,«;||2 + A^+i,,(l-a)||u7|p. (6.33) 

We determine a by writing 

^N+i,z{^ — a) = Xi^z , 
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hence 

a = 1 - . (6.34) 

So we have 

HQ ||2^ 2 TV ^ „^N+l,z N m otx\ 

\\dzw\\ < -mB,n <2— rriB^. (6.35) 

a dN,z 

In the semi-classical regime where we are, this leads to the existence of a 
constant C such that 

\\d,w\\'^ <Cm^^^. (6.36) 
Using in addition the diamagnetic inequality, we obtain 

\Mw\\\l<Cml^. (6.37) 

As in the other case, we obtain from Sobolev's Inequality the control of 
w in norm 



Me < Cim^,n)Hl + < C^Kc)^ , (6-38) 



where we have used that 6^ >> 1 in the semi-classical regime. 
Using Holder's inequality, we obtain 

lkll4<C(m^,f,)^(0-^- (6.39) 
We now have all the estimates needed to mimic the proof of case A. 

We start from 

^(^) > Ai,^ + ^ij(njv^) - 4c/ y in^v^l VI dxdydz . (6.40) 

We have now to control the third term in (16.401) by the second term. This is 
done like in case A in the following way : 

Ag J \UN^f\w\ dxdydz < 4^| lU^v^l |3 | |w| I4 
We now use 

SB{IlN-^)>m^^^\\UN-^\\t, (6.42) 
which together with fl6.29p leads to 

N 

m% < C{1 + ^)Sb{II^^) . (6.43) 

z 
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This leads to 

Ag j \Ii^^^\w\ dxdydz <C2g-^{ml^^)"^{5^)-l{l+'^)^^^ (6.44) 

Using this control, fl6.29p . fl6.4UI) and fl6.42p . we have obtained the detailed 
proof of fl6.17p in the general case. 



6.4 On the minimizers of 8b' 

The next proposition is rather standard and refers to the case = 1 but has 
its own interest. As a corollary, this will yield an upperbound for mB,n- 

Proposition 6.10. 

The minimizer of £b over the normalized ip 's is unique (up to a multiplicative 
constant of modulus 1 ) and radial. 

Proof : 

We first observe that if ^/^ is a minimizer then is also a minimizer. Con- 
sequently, we start considering a non negative minimizer. 
Now \ is solution of the corresponding Euler equation and by the Maximum 
Principle, cannot have a local minimum. Hence ip cannot vanish and we 
can write 

Comparing Sb{iP) and Sb{\iP\) we get 

|V^| = |V|V'| I a.e. 
and this implies that a is constant. 

So we can now assume that ip is a real positive minimizer. The Euler- 
Lagrange equation reads 

- Alp + uj_i_^r^4j + gl^pl'^i; = (6.45) 
for some real Lagrange multiplier A. 

Let another positive solution (with ||0||l2 = 1) of the Euler-Lagrange 
equation for a possibly different G M : 

- A0 + a;±V0 + ^|0|V = (6.46) 
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Possibly exchanging the roles of ip and 0, we can w.l.o.g assume that 

A > /i. 

Let us consider, for some a > to be determined, the rescaling 

y) = \fa u{\fax, \fay) , 

we get for u the equation 

—Am H — r M + (yf m m = — m . 

Oi 

We now choose « = ^ which leads to 

- l\u^u)^ (^^r'^u^ g\uS^u = \u. (6.47) 
We can now compare u and ^l). Let us introduce v = ^ which is a solution of 

-div (V-^Vt;) + gvip^{\v\^ - 1) = cJxV^^t;^/'^/^ _ _ 
Multiplying by (f — 1)+ and integrating we obtain 



A^ 

^'^\V{v-l)+\'^ + gv^'^{v-l)\{v + l) dxdy = j uj±'^r^v^^{l- —) dxdy . 



With our assumption on (A, /i), this implies the vanishing of {v — 1)+ almost 
everywhere hence f < 1. 

This can be reinterpreted as m < ^ and the normalization implies u = ip. 
The lemma is proved. 

Finally, one can construct a radial positive solution (by minimizing Sb 
over the radial functions). This gives a solution (f) of the Euler equation which 
is also strictly positive. ■ 

Observing that, if ip is radial, we have that £b,q.{.4^) = £b{4^), this propo- 
sition has the following interesting corollary. 

Corollary 6.11. 

We always have 

inf Sb,^ ■■= mB,n < mB ■ (6.48) 



^In full rigor, we should consider a sequence Xn{y — 1)+ where (xn)n is a suitable 
sequence of cut-off functions and take the limit (See |BrOs| ). 
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6.5 Lower bounds in the TF case (A^ > 1) 

We start from a minimizer {ipe,±)i- Due to the normalization, there exists at 
least one j such that 

ll^,,±ll> ^ 



Then we write (neglecting the kinetic part) 



NT /N~l 



2 "'•x.y \j=0 



When expanding (j:^ i,^ (.M,,,(..y))' , the mixed tenn. a.e exponen- 

tially small (see Subsection 17. ip in comparison to \ \ipj,±\ |^4, hence we get, 
for some a > 0, 

NT 

<,n > li^'-^') I A^,A'+9 j_'^i^o{zfdz{ j {^P,,^ydxdy) (1-exp-^) 

We now use (17.41) . to obtain 

<,n > U^'-n')Jr'\^,M' + 9I^M^rdz{J^Pl^dxdy){l-exp-^) 
= liu^ - J^2|^^. j2 + g(^j;^^^dxdy){l - exp -f ) 

> (1(^2 _ ^2) J ^2|^^. + ^l^dxdy)) (1 - exp -f ) 

> ^(1 - exp -f ) inf^,||^||=i (l(a;2 - Q^) J r^j^p + g{J rdxdy)) . 

One can then use the asymptotics obtained in the proof of (I6.16P to get, 
under Assumption (I6.15p . the existence of Cn,i3 > such that, as e tends to 
and g to oo, 

m%>y^oo^y^. (6.49) 



7 Tunneling effects for the non-linear models 

This is only in this section that we will exhibit the role of these localized 
(A^T) -periodic Wannier functions. 
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7.1 Towards the DNLS model. 



7.1.1 Preliminaries 

We have already proved for any > 1 the rough estimates on allowing 
to justify the longitudinal model. We have established or announced better 
asymptotics at the price of stronger assumptions. 

Our aim in this section is to discuss possible asymptotics for in the case 
when > 1, which will involve the tunneling effect. Although we have no 
final result on this part, we would like to prove how we reach a familiar model 
considered by the physicists : the DNLS model. 

In particular we will describe in Proposition 17.71 under which assumptions 
one can get a simplified model. 

We consider on ] — ^[ the (A^T)-periodic problem for the operator 
— ^ + We{z) with We{z) = w Satisfying Assumption (11.11) . Here we 
always work in the semi-classical regime. 

The starting point in this subsection is that we replace the issue of min- 
imizing £^'^'^ on the (A^T)-periodic L^-normalized functions by restricting 
the approximation to the eigenspace ImvTTv associated with the first A^ eigen- 
values of the linear problem. 

7.1.2 Projecting on the eigenspace ImTTjv 
Our aim is to analyze the reduced functional 

j=0 

where £^'^'^ is in fact £^ given in fll.35p with the explicit notation of the 
dependence of the parameters and the ip^ are the (A^T)-periodic Wannier 
functions. When A^ = 1, the error which is done has been estimated in 
f lS.lSp under the assumption that ge^^ is small, i.e. (11.491) . Replacing in the 
argument the projection on the first eigenspace by vtat, the same result holds 
true for A^ > 1. 

We now concentrate our discussion to the model obtained after this first 
approximation. More specifically we are interested in the asymptotics of the 
infimum of this functional. 
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Here natural approximations of this reduced functional appear. Each of 
these approximations gives a corresponding approximation of the infimum of 
the reduced functional, which is defined by : 

m^/'\e,g):= M £^-'^'-\c) . (7.2) 
{c| Ef="o' k,P=i} 

Let = Xi^z be the bottom of the (A^T)-periodic spectrum of on 
[~^^^] (with minima). So strictly speaking, we can start the anal- 
ysis of this first approximate model only under Condition (11.491) . 

Proposition 7.1. 

Under condition (11.491) 



m^{e,g) = m^'^'\e,g) + 0{gle-'^) . (7.3) 

One can nevertheless imagine that the information obtained in the next 
subsubsection is valid in a more general context (maybe by choosing other 
localized Wannier functions). We now analyze various approximations of 

mA^°\e,g). 

7.1.3 Neglecting the tunneling 

Neglecting the tunneling effect, we are lead to the minimum of the functional 

(N-l \ N-l „iVT 



j=Q J j=0 



over the c's such that 



N-l 
j=0 



Observing (see |DiSj| ), that 

NT T 

[ " \i^^{z)\'dz= r 0,(^)4 rf^ + a(exp-f), (7.4) 

J_NT J— I. 



2 2 



where 0i is the groundstate of the T-periodic problem, the minimum of this 
approximate functional, which is attained for Cj = N~2^ is 
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mT^ = \i,. + i \i^^{z)\Uz. (7.5) 



NT 
2 



So as a first approximation, we have obtained 
Proposition 7.2. 



mf^'\e,g) = \,., + ^( / ' {z)\' dz) + (?+ 1) a(exp--) 
iV J_NT e 



or 



(e, g) = + M^V dz) +g 0(exp + 0(exp --) . (7.6) 

Tlie definition of O is given in fll.28p . If we apply tfiis result to our 
context witli ^ = ujj^g, tliis yields information on the behavior of m^'^°'* 
independently of Assumption (ll.49p . 

7.1.4 Taking into account the tunneling 

If we keep the main tunneling term, we get the following more accurate 
approximating functional 

■■= Ai (Ef=o' Ic,f ) - (Ef=o' 9Tt) + ?(E;=o' Ic.f) (/i \i'^{^)\'dz) 

(7.7) 

Here r is the hopping amplitude introduced around (13.121) . Ai is the low- 
est eigenvalue corresponding to the Floquet condition /c = y for the linear 
problem on ] — |-[, which is exponentially closed to Ai and we take the 
convention that Cat = cq. 

The quadratic form corresponds to the approximation in the first band : 



'N-l \ /N-l 

3c^Ai (5^|c,f -r3fj K^c.c^) (7.^ 

3=0 / \ j=0 



which can be shown to be correct modulo 0{exp — 



25 ^ 



Remark 7.3. 

This time the minimizer could depend on 'g !! This is the kind of problem 
which is analyzed in IKMPSf and in Subsection \7.3 
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Discussion about the justification of g^'^'^'^'^^ 

One can wonder why we forget some terms in the computation. Let us 
do this more carefully. To be consistent with what we forget in the lin- 
ear case (terms of order (9(r^)), we show first that one can approximate 



{lJ^\E^oC,^n^)\'dz) by 



2 
NT 



/ NT 

(/_klEf=o'c.<(;^)N^ 

, NT 



2 

/ NT 

+ p=o (dc.f + l9+iP)(c,C(-TT) + c,+iC(-y) {j:^^^{z)\^^{z)\'-^f{z)dz) 

(7.9) 

This first approximation is based on the following lemma. 
Lemma 7.4. 

tP^{zf^^{zfdz = 0(exp - — ) . 

NT 6 
2 

This is based on the property that, for all rj > 0, there exists C.^ such 
that 

\^^{z)\ < C.exp I exp --^dT/iz) , (7.10) 
where d'Xfi z) is an even function such that 

dTg\z, 0) = 2 / ^^Mf)dt , for ;2 G [0, T[ , 
Jo 

and such that d'^g'^{z, 0) is increasing for z > 0. 

On the contrary, this is a priori uncleail§ why one could forget terms like 

r = 9 I ^^{zf^nz)dz. (7.11) 

/ NT 

(where we recall that w is even by Assumption (11.11) and that this implies 
-00^ even and real). This term is a priori of the same order as r. We have 
indeed 



®We use here the assumption that the potential and hence ipQ is even. We recah also 
that the tpj are real. 

^In |KMPS| . p. 5, between formulas (18) and (19), the term r is discussed; see also p. 6 
around formula (20). 
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Lemma 7.5. 



NT 



I " ^^{zf^n^) dz = 0(exp --) . (7.12) 
/ NT e 

Due to the decay estimates (17.101) for these {NT)- Wannier functions, the 
term to integrate in (17.121) decays hke 

d[e^V-\{U7f{z)+dAf{^-T)) 

so the main contribution comes from the origin and has the same size as 
exp — ^ . 

So it is necessary to be carefuf^. if one wants to neglect r. 



Let us now try to estimate J {zY'ipi{z) (iz as e — more precisely. 

Heuristically, one can try to use a WKB approximation, this is available for 
TpQ in the neighborhood of but unfortunately, we do not have a good WKB 
approximation of tpi {z) close to the origin, as observed in Subsection 13.31 
(see (I3.23P ). So we have no obvious main term for the asymptotic behavior 

NT 

of j_NT_ i'o {z)^ip^{z) dz. A reasonable guess (which is implicitly used by 

2 

the physicists) should be to suggest the following conjecture. 



Conjecture 7.6. 

r = ?o(r), (7.13) 

as e ^ 0. 

The weaker mathematical result, which is obtained from Lemma (7.51 by 
the considerations above using Helffer-Sjostrand techniques |DiSj| , is the fol- 
lowing proposition. 

Proposition 7.7. 

Under the assumption that there exists rj > such that, 

0<gexp-<l, (7.14) 
e 



10 



We thank M. Snoek for kindly answering our questions on this problem. 
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then 



m 



7V,(0) 



m 



TV, (2) 



+ r 



(7.15) 



holds. 



This gives a motivation for the analysis of the DNLS model of |STKBj 
(with an extra term in A ^^=q l^jP)- 

If we consider the (iVT)-periodic Floquet problem, we arrive naturally to 
questions analyzed in |KMPSj (16-17-18), and the remark after (21) in this 
paper. 



7.2 On approximate models in case B : towards Snoek's 
model 

Using the basis of the (A^T)-Wannier approach), we can consider 8^ intro- 
duced in (11.421) and consider the decomposition 

N-l 

(^/'CX, ■ ■ ■ ,'0iv-i,±) :=^b'(^o,±,--- ,'tpN-i,±)+g\\ {z)i)j^^{x,y)\\[i . 

j=0 



We now use various approximations related to the analysis of the z- 
problem ((A^T)-Wannier functions). We get 



'(^o,±, ■ ■ ■ ,^iV-l,±) 

~ s J2fJo^ I \i^j,±\? + 1 Ef=o^ + (V'i.x, V'i-i,±)) , 

and 

So the approximate functional becomes 



^N,appro.^^^^^^^^^ = ^fj; f^, ilm^J' + r(x,y)|V^,,^(x,l/)H dxdy 

+sEpo^\\i^,,±\\' 
+9J2f=o Il^j,±lli4, 

(7.16) 
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which should be minimized over the {ipj^±)j such that 

Af-l 
j=0 

This is the model described by Snoek [Snj . 

Starting from this model, one can, depending on the size of the various 
parameters, come back in some case to the situation when {4'j,±)j is of the 
form Cjijj±, with J2j=o = 1- In this case, we come back to the results of 
the previous subsection. In the other cases, the problem seems completely 
open. 

7.3 Spatial period-doubling in Bose-Einstein conden- 
sates in an optical lattice 

Here we mainly follow [MNPSj . We look at the discrete model for which 
these authors refer to [STKBj . 

The Hamiltonian for the discrete model is formally definecJ^ as 

j j 
where c = (cj) G ^^(Z; C), r > 0, / > 0. 

This is a particular case of the so-called DNLS modelll (Discrete Non 
Linear Schrodinger model). But we will immediately reduce our analysis to 
the (A^T)-periodic problem and Floquet variants of this problem. 
We will restrict the sum above to j = 1, . . . , A^, where is a fixed positive 
integer. But for defining the tunneling term (cjcj^i + CjCj^), we use the 
Floquet condition 

cn+i = exp{ikN) ci . (7-17) 

Remarks 7.8. 

• IMNPSf looks at the particular case N = 2^ . 
""^^after substraction of a term in the form fX^j 

-•^^In the general case one should add the term ej|cjP when the trapping is present 
the authors propose Cj = ujj'^ with oj proportional to uj^. 
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• When k = 0, this is the natural {NT) -periodic problem. 

• Note that it is the standard Floquet problem only for N = 1. 

• Note that we forget the term A Ylj which corresponds only to a shift 
of the energy. 

Hence we would like to minimize H^'^{c) 

N N 

H'^'^ic) := -T J2ic]c,+, + c.cJT,) + I J2 Ic.f , (7.18) 

over the c's normalized in by 

l|c|r = iVc, (7.19) 
and the fc-Floquet condition f l7.17p . 

Moreover for future use we introduce the strictly positive parameter 

u = NJN > . (7.20) 

Remark 7.9. 

In the preceding sections we were taking Nc = 1. Up to a change of the 
parameter r, one can always reduce to the general case to this situation. A 
difference could occur if we take u fixed and N +oo. 

We will then be interested in the analysis of the energy per particle 

E{T,I,u,N,k) = ^ inf H'^'^ic). (7.21) 

-''C 11*^11 — -"c 

Writing 

Cj = exp i kj Qj (7.22) 

then from fl7.17p 

9i+N = 9i ■ (7.23) 

The case when = 1 corresponds to the usual Floquet condition. 

Writing that c is a critical value of H^'^, we get that for some (which is 

called the chemical potential or the Lagrange multiplier in mathematics) 

2I\cj\^Cj - T (cj+1 + Cj_i) = fiCj , (7.24) 

for j = 1, . . . , A^, with condition fl7.17p . This becomes in terms of the gj^s 

'^I\9j?gj - (expik gj+i + exp{-ik) gj^i) = figj , (7.25) 

for j = 1, . . . , A^, with the A^-periodic convention fl7.23p . 
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The case N = 1 

In this case, we have simply one equation 



2/|^i|2^i-2tcos(/c)^i = //^i, (7.26) 

with 

= (7.27) 

We find immediately that 

II = -2r cos(A;) + 2Ip (7.28) 

and 

gi — i^expi^i . (7.29) 
The energy per particule E is then equal to 

E = -^H^'^c) = -2t cos(A;) + lu . (7.30) 

For this choice of N, we have recovered exactly what we have found in the lin- 
ear case 7 = 0. The effect of the non-hnear term just creates a A;-independent 
shift of the energy. Note also that as a function of k, the energy is minimal 
for A; = 0, which is the periodic case. 



The case iV = 2 

In this case, we get, using the periodicity assumption, the following system 
of equations 

2/|«/i|Vi - 2tcos(A;) ^2 = ligi , 
-2t cos(A;) Qi + 27|g'2p5'2 = fJ'92 , 

with the normalization condition 



(7.31) 



\gi\^ + \g2\^ = N, = 2u . (7.32) 

We write 

9j = \9j\ expi(/7j- , for j = 1, 2 . (7.33) 
A suitable combination of the two lines gives 

2/ (Ig^il^ - 1^2^) = 2TC0s(/i;) ( |^expi((/?2 - V'l) - exp -i((^2 - V'l) ) ■ 

\92\ / 
(7.34) 
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If we observe that the left hand side is real, we get 

= 2rcos(A;) ( ^ + |^ ) sinfy^a - V'l) • (7.35) 
The real part of (17.341) is in any case given by 

I - \92\^) = Tcos{k) (|5(2p - Is-ip) - — I — rCos(v32 - <^i) . (7.36) 

We meet three cases. 
Case 1 

We first observe that the solutions corresponding to the case = 1 are 
recovered by taking ipi = and \gi\ = \g2\. 

Another family of solutions is obtained by taking ipi = y^a + vr and \gi\ = \g2\. 
This corresponds to an "antiperiodic" solution over two periods. 

The solutions such that \gi\ = \g2\ = seem to be simply deformations 
of the case 1 = 0. 

At least for I small, this is indeed also a consequence of the implicit function 
theorem when cos k ^ 0. 

The energy per particle E is the same as for A^ = 1. 

Case 2 

If we assume that 

\9i\7^\92\, (7.37) 

and 

r cos A; 7^0, (7.38) 
then the previous necessary conditions become first 

sin(v92 - v^i) = 0, 

hence 

ip2 = (pi mod{Zn) , 
and secondly (using the first one) 

2/|^i||^2| = ±2rcosA;. (7.39) 

Using the normalization of gi and g2, we get as a necessary condition 

|rcosA;| < . (7.40) 
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If these conditions are satisfied, we can find in function of the sigio of 
cosfc, a unique pair gi, g2 (up to a multiplication by e*^). Coming back to 
the initial system of equations leads to the determination of fi which is given 

by 

/i = 4/i/, (7.41) 
and of the energy E which is given by : 

E = j^cos^k + 2Iv. (7.42) 

Case 3 

It remains to consider the degenerate situation. If 

cosfc = 0, (7.43) 

i.e. if 

\k\ = ^ 



\k\ = l, (7.44) 



Then 

91 = 1^ exp iipi ,92 = 1" exp iip2 , (7-45) 
is a solution for any pair {ipi, ip2)- The corresponding /x is given by 

fx = 2Iu, (7.46) 

and 

E = Iu. (7.47) 

This is the same energy that the one found for the usual Bloch state (with 
the same k) but note that the gj are no more 1-periodic {gj+i 7^ gj). 

Question 7.10. 

It is unclear in the discussion what is the status of k. Are we interested in 
minimizing over k ? But in this case k = and N = 1 seems optimal in the 
sense that they give the lowest E. 

Remark 7.11. 

An interesting problem, which is discussed in IMNPSf . is the analysis of the 
stability, looking at the linearized corresponding problem. 

Remark 7.12. 

Of course the analysis of more general N 's would be quite interesting. A few 
numerical results are given in JMNPS^ corresponding to N = 4. 



13 



It seems that in our problem we have / > and r > 0. 
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8 Other optical lattices functionals 

In this section, we discuss the choice of analyzing periodic boundary condi- 
tions in the z direction and the possibihty of stating the problem differently. 
We compare the (A^T) -periodic problem to the T-periodic problem and dis- 
cuss shortly the question N -|-oo. 



8.1 Summary of the linear case 

We summarize what we have obtained in the linear situation. Different tech- 
niques can be used for determining the bottom of the spectrum of and 
then of if, but the ground state energies always coincide. 

(i) Minimize the functional 



over L^-normalized ■i/^'s in C^(M) (or in if^(R)). In this case, the 
minimization gives the ground state energy but there is no minimizer 
in the form domain of the operator ! 

(ii) Minimize the functional 



over the I? normalized C°° T-periodic functions '?/''s (or on H^'^*^^) . 
Here we integrate over one period ! The minimization will give the 
ground state energy of the periodic operator and the minimizer of the 
functional will be the ground state. 

(iii) Minimize the functional 



over the L^-normalized C°° (A'"T')-periodic functions (or on H^'^'P^'^). 
Here wc integrate over periods ! The minimization will give the 
ground state energy of the periodic operator and the minimizer of the 
functional will be again the T-periodic ground state. 



(8.1) 





2 



(8.3) 
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(iv) Minimize over k G [0,27r/T[, the infimum of the functional 

= (^l^'(^) - + W,{z)mz)\'^ dz (8.4) 

over the L^-normahzed T-periodic functions (or on H^'^'^^'). Here 
we integrate over one period ! The minimization will give the ground 
state energy of the periodic operator (i.e. = 0) and the minimizer 
of the functional will be the periodic ground state (corresponding to 
k = 0). 

(v) Minimize over the space generated by the Wannier functions (identified 
with ^^(r, C)). This corresponds to the reduction to the first band and 
leads to the analysis of 

£2(r,C) 3c = (Q)ter ^ J] Ai(£-m)Qc;;;, (8.5) 
with Ylii = 1- 

(vi) Minimize over the space generated by the (iVT)-Wannier functions 
(identified with i'^(T^ , C) = C^). This corresponds to the reduction to 
the spectral space attached to the first eigenf unctions of the (NT)- 
periodic problem living in the first band. 

We emphasize, that, in each case, we get the same ground state 
energy. 

As already mentioned, the 3-dimensional linear case introduced in (11.161) is 
easily reduced to the one- dimensional case Hz- 

Remark 8.1. 

The reader could be astonished that we discuss only the case = 0. This 
is a current assumption in the physical literature. Mathematically, there is 
a dramatic change in the nature of the spectrum. The problem could become 
quite difficult in some regimes of oj^ and this will not discussed in this paper. 
Let us nethertheless make a few comments. When uo^ 7^ 0, the spectrum of 
Hz on the line becomes indeed discrete. By monotonicity, the bottom of the 
spectrum is above Xi^z- One can also get an upper bound by computing the 
energy of a suitable quasimode (or more simply of the ground state of the 
linear problem) but this can only be good in some asymptotical regime. 
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8.2 The 3L>- functionals 



Let us consider the fully non-linear problems and try to implement some 
of the results obtained for the linear models. We take Uz = 0, we have to 
analyze (see (11 .ip ) the functional 

^^gn(^):= 

J^^ Q|(V,,,,, - zfi X r)vl/(r)|2 + (\/(r) + W,{z) - ^n''r^)\^iv)\'' + (^[^(r)^^ dxdydz 

(8 

We denote by the natural maximal form domain of the form (5n(^) in 
L2(R3)^ that is 

Vn = {ue H\R^) , xu G ^^(R^) , yu G L'^{R^)} . (8.7) 
We denote the intersection of the L^(M^)-unit ball with Vq : 

Sn = {^ eVn, 11^1^2 = 1}. 
We call the infimum of this functional 

En := inf QnW . (8.8) 

Because there is no harmonic trapping in the 2;-variable the choice of the 
condition of normalization in L^(]R^) is questionable. 



The T-periodicity of the optical lattice in the z-variable suggests to con- 
sider other functionals, where we integrate over K^^^x] — ^,^[ (or more 
intrinsically over M^. ,^ x (M/TZ) ), or over ffi^.^yX] — ^[ for some integer 
N, and where the variational space has to be defined suitably (periodic con- 
ditions or Floquet conditions). We refer to Subsection l8.lt for the discussion 
done in the linear case. In the non-linear case, this has led to the introduc- 
tion of the "periodic" Bose-Einstein functional (see f ll.lSp ) We will denote 
by Qn^'^ the functional obtained by integration over periods (see f ll.lOp ). 
We call T>^^, the natural maximal form domain of the form ■ It corre- 
sponds to the distributions in Hl^^(M.^), satisfying fll.l3p and such that, 
the restriction \E'r to K^^^x] — -j, |-[, satisfies : 
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We note that this functional has clearly a minimizer in Sf^^ , where 
cSr = {^G^^r, / \^{x,y,z)\'dxdydz = l}. (8.10) 



We denote this infimum by 



In the spirit of the Floquet theory, one can also be interested in the analysis 
of the Floquet Bose-Einstein family of functionals, defined for G M by : 

-/«i.x,-X.[(il(V.,.-^M^P-i^V|vl.p 

+ ||(a, + ik)^\^ + {V{r) + W,{z))\^\'') dxdydz ^ ' 

+g ki^^y^yTT^m^ dxdydz , 

where \l/ satisfies 

^(x,y,z + T) = ^(x,y,z). (8.13) 

We call T)^'^'^^ the natural maximal form domain of Q^"^^ which is actually 
independent of k and VL : 

^Floq 'pfPer 

^n,k ~ ■ 

We note that we have here a family over k of functionals. Each of these 
functionals has a minimizer in S^^""^, where 

=: {M/ G Pgr , / \^\' dxdydz = 1}. (8.14) 

It is natural to be interested in the quantity 



with 



Eg°'':=infES^ (8.15) 



<r:= inf gSr(*)- (8.16) 



f e<s, 



Floq 

n,k 



Remark 8.2. 

In the physics literature, this corresponds to a ground state of a condensate 
at rest in the frame of the optical lattice. The energies -E^'^ describe states 
of the system where all the atoms move with respect to the optical potential, 
at fixed velocity, giving rise to a constant current equal to k. Experimentally, 
this is achieved by moving the optical lattice with respect to the condensate. 
We refer to jKMPS^ and the references therein. 
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We would like to compare the various functionals. 

Because g > 0, all the functional are bounded from below, admit an 
infimum and a lower bound is given by the analysis of the linear problem 
corresponding to g = 0. But there is probably no existence of a minimizer 
for the Bose Einstein functional. So we will have to consider minimizing 
sequences or approximate minimizers using possibly other models having 
minimizers. We refer for this to Subsection 18. 1[ The new point is now that 
it is unclear that the various ground state energies obtained by the different 
procedures coincide ! 

Finally, it could be interesting to compare directly the infimum of the 
three functionals above and then to find good approximations of these infima. 

At the moment we have shown in the "linear" section that : 

En{g = 0) = ^^(^7 = 0) = E^^'^^g = 0) = mfa{Hn) = \. + oJi. . (8.17) 

We have also mentioned, observing the monotonicity of the functionals with 
respect to that, for > 0, 

En{g = Q)<En{g). 

E'-{g = 0)<E^-\g), (8.18) 
and E^^"\g = 0) < E^^'^^g) . 

This is further developed in Proposition 18. 3[ 



8.3 Comparison between the Floquet Bose-Einstein func- 
tionals and the periodic Bose-Einstein functional 

The argument which follows is only correct when i7 = 0. The idea is that 
using the Kato inequality, one has 

Q"B^(l^l)<Qr''W, VA; (8.19) 

In the periodic case, the minimization over complex functions leads to 
the same infimum as in the real case. 
We then obtain 

hence 

^Floq^-^per_ (8.21) 

This seems to suggest that, when = 0, there is no interest to introduce Flo- 
quet conditions, if one is only interested in the determination of the ground 
state energy. 
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8.4 Comparison between the periodic Bose-Einstein 
functional and the Bose-Einstein functional on 

Here we again work for general Q's satisfying fll.Sp . We can naturally try 
what was working in the linear case. If \E'^^'"(x, z) = 'il)per{x^y)ipi{z) is the 
periodic minimizer of the linear problem, i.e. satisfying 

Hamper ^ j^Un^per ^ 1 1 ^^^11 ^.(Mi,, x]- f ,f [) = ^ 

and 

W^'ix, y,z + T) = y, z) , (8.22) 

with i?'*" = uj± + )^i,zi we can use as trial function for the Bose Einstein 
functional 

^R:=CRXR{z)^'"'{x,y,z), (8.23) 
where cr > is determined by the condition 

4/ \xR{z)^!^''{x,y,z)\Uxdydz = l, (8.24) 

Here i? > 1 is a free parameter which will tend to +00 and xr is a function 
with support in [—RT — T, RT + T] equal to 1 on [—RT, +RT] whose first 
derivative is independent of i? > 1. 

It is immediate to see, from the normalization chosen for vI/P'^'^^ that 

Cr ~ — , as i? ^ +00 . (8.25) 
v 2-R 



We note that these trial functions are radial in the (x, y) variable and are 
independent of fl. 

The main point is to observe that we have 
/]g3 \'i!R{x,y,z)\^ dxdydz ~ 2i?c|j /j.2^j_t_tj |**''""(x, z)|*^ dxdydz 

~ 27? /R2x]-f ,f [ y, z)\^ dxdydz . 

(8.26) 

Hence 

lim / \^R{x,y,z)\'^ dxdydz = 0, (8.27) 
R^+'^ Jr3 

We also obtain easily that, as -R — >■ +00, 
lim / \{V-inxr)^R{x,y,z)\'^ = [ \{V-znxr)'ifP^' {x,y, z)\'^ dxdydz 



(8.28) 
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and 



lim / (l/(x,7/,z)-^nV)|*^(x,?/,, 

= [ {V{x,y,z)-ln\')\^'^^'{x,y,z)\'dxdydz. (8.29) 

So we obtain that 

lim Qa{^R) = Ef,'^{g = 0). (8.30) 

it— *+0O 

One can also observe that everything is actually ^2 independent. 



Combining with what we have verified for the linear case, we obtain that, 
for g > 0, 

El^^g = 0) = E{g = 0) < E^{g) < lim Qa{^n) = E^^^ig = 0) . (8.31) 

i?— »+oo 

So we have proved the 

Proposition 8.3. 

For all g >0, allO <9, < u±, 

E{g) = Ea{g) = E'^''{g = 0) = Ea{g = 0) = E{g = 0) . (8.32) 

The conclusion is that if we look at the Bose-Einstein functional on 
the infimum of the functional restricted to L^-normalized states is indepen- 
dent of > and fl and is immediately obtained by the ground state energy 
of the Hamiltonian attached to the case g — and = 0. 



The result is of course also valid for the 1-dimensional problem and is 
independent of any asymptotic analysis. 

Remark 8.4. 

Another natural physical problem would be to analyze the quantity 



liminf — 



( \ 

inf Q^-'^{^) 

\!^^\^\''dx=N, j 
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or 



lim sup — 



( 



\ 



inf 



I 



J ^ |*|2dx=JV, 

where we compute the energy by integrating over N periods and where 

NJN = V 



(v fixed). Upper bounds for this model are the periodic models with g replaced 
by gv. We met already this problem in Subsection \ 7. for discrete models. 



8.5 Comparison between the (A^r)-periodic problem 
and the T-problem 

In this subsection, we pursue the analysis of the hnks between the (NT)- 
periodic problem (A^ > 1) and the T-periodic problem. We recall from Sub- 
subsection [2l2l3] that, for the (iVT)-periodic problem in ID, the ground state 
energy is Ai^^. Moreover, in the semi-classical regime we have a packet of N 
eigenvalues which are exponentially close separated from the (A^+l)-th eigen- 
value by a splitting which satisfies (11.271) . In case A, a natural question is : 



Question 8.5. 

Is the minimizer of T-periodic as in the linear case ? 



When the answer is yes, we immediately obtain that 

"^A(e>^) = "^A(e, ^) (8.33) 

and we can directly use what we have done for proving Theorem 11.31 by 
replacing 5^ by f . 

To our knowledge the answer to this question is unknown, so it is natural to 
look at simpler models. 

In case B a natural question could be 
Question 8.6. 

Is the corresponding minimizer in this reduced space T-periodic ? 
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When the answer to this question is yes, we have seen in (16 ■4p together 
with the discussion around (I5.43P (with the additional assumption that the 
minimizer is periodic) that 

We could then use what we have used for the proof of Theorem 11.51 This 
should work in the case when g and VL are small. 

Actually, we could have asked more directly the following question 
Question 8.7. 

Under which condition on g and Q is the minimizer of the [NT) -periodic 
initial problem T -periodic ? 

If it is the case, we get immediately 

Er^'ig) = Ef,-ij^) (8.34) 

So we can directly relate the treatment of the (A^T)-periodic problem to our 
preceding studies, without any use of (A^T)-periodic Wannier functions. 

The general answer is unknown. One suspects by bifurcation arguments 
that it is true for g and Q small enough, but the physicists seem to wait for 
the other situation. This should in particular be the case for sufficiently large 
rotation Q (see for example |CorR-DY] ) or in the case with (A^T)-Floquet 
conditions (see the discussion in Subsection 17.31 on an approximating model). 



A Floquet theory 



We follow the presentation of |DiSj| (p. 160-161), who are actually dealing 



with a more complicate situation. If we take as lattice T = ZT, the dual 
lattice r* is T* = ^Z. For u G S{R) and e M, we put 

Uu{y] k) = exp Z7 k u{y — 7) . (A.l) 
7er 

We notice that L{u{y; k) only depends on k modulo the dual lattice V* so 
Uu(?/; k) is well defined on M x (R/F*). 
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For k e M/r*, we put 

V'^, = {ueV'{R); u(y + 7) =expi7fc u{y)} , (A.2) 

and 



nk = {ueLlM)r^v',■, [ 

Je 



uiy)\' dy < , (A.3) 



where £" is a fundamental domain of T (for example E = [—^, ^[). 

When k = 0, Hq denotes simply the space of T-periodic functions in Lf^^{R). 

From ( lA.ip we see that 

Uui-;k)enk, (A.4) 

and if we view Tik as a bundle over M/F*, we can view Uu as a section of 
this bundle. We write Uu G C°°(M/F*; Hk). 

Now, if u G L^(M/F*), one can expand it in a Fourier series : 

v{k) = ^^^(7) expijk, (A. 5) 

where 

^(7) = — / 6xp — ^7 k v{k) dk . (A. 6) 
2vr Jo 

We have the Parseval formula 

f r \v{k)\^dk = Y^\v{^)\\ (A.7) 



Now, for any y G M, we can view ( 1A.1I) as the Fourier expansion of ■). 
so ( 1A.7P gives 

|W«(y;A;)prffc = 5^Ky-7)p. (A.8) 
Integrating over y G [— -|^], we get 



^ /I |Wu(y;A;)pdA;j dy = ^J^l\u{y-j)\' dy = £^ \u{y)\' dy . 

(A.9) 



2 ' 2. 

+00 



27r /_T „ ^ 

2 V"^ " / 7er 2 
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So U can be extended to an isometry from L^(R) into L^(M/r*; Tik)-, which, 
as an Hilbert space, can be also described as the space of the w's in Lf^^{M?) 
such that 

with the norm, whose square appears in the left hand side of ( jA.91) . i.e. 



\v{y]k)\'^ dydk. (A.ll) 

[-|,|:]x[0,2VT] 

Like in the standard analysis of the Fourier transform, we have now to ana- 
lyze the surjectivity property and the construction of an inverse. 
If {y-ik) I— > v{y]k) belongs to C°°{R/T*;T-Ck), we can write the Fourier ex- 
pansion of V with respect to the second variable : 

viy;k) = ^v^{y) expi'yk, (A. 12) 

with 

27r 

v-y{y) = — / exp — 27 v{y;k) dk . (A. 13) 







271 

Using the first line of ( lA.lOp . we see that 



%{y)=vo{y-'-r)- (A. 14) 

Hence, when v is smooth, we have 

v = l{yVv, (A.15) 

where 



^v{y) = — / v{y; k) dk = Vo{y) ■ (A.16) 







Using (IXTHll . dsn and (lAl2l) . we obtain 



l|w^lli.(M) =T.,^T!^T\My-i)?dy 

= i:,^vpT\v,{y)?dy (A.17) 

2 ^ 
= ll^llL2{R/r*,Wfc) • 

This shows that W is an isometry of /.^(IR/r*; 7^^) into L'^{M). Hence W is 
a bounded right inverse for U. Since U is an isometry, we conclude that U is 
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unitary and that its inverse is W 

dy 



Now, when considering P := —4^ + W{y) with W T-periodic, we see that 



UPU-^ = Pkdk, (A. 18) 

J[0,2tt/T] 

where, by definition the right hand side in (1A.18P denotes the selfadjoint 
operator Q on L'^{R/T*;Hk), with domain L'^{R/T*;Hl), which is given by 

Qv{y;k)= (^PM-,k)){y). (A.19) 

Here above 

nl = {uenk, m(^) G Hk for \i\ < 2} , 

and Pk is the selfadjoint operator on Hk associated to the differential operator 
-£2+W{y) with domain Hi- 

Finally, we note that for k G M/F* , Pk is unitary equivalent to the operator 

Pk = exp—iyk o P^ o expiyk, (A. 20) 

which is now acting on Tig (corresponding to the T-periodic fuunctions in 
Lf^^). More explicitly it takes the form 

d 
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